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НАПРЯЖЕННО-ДЕФОРМИРОВАННОЕ СОСТОЯНИЕ 
ЗАРЯДОВ РДТТ В ОБЛАСТИ ДОННОГО ТОРЦА 

В. А. ДУНАЕВ, А. А. КАШИРКИН, О.А. ЕВЛАНОВА 

ГНПП «Сплав», Тула, Россия 


АННОТАЦИЯ, Разработана комплексная математическая модель, содержащая мо- 
дели одномерного квазистационарного течения, двумерного течения газов в области 
торцов заряда и напряженно-деформированное состояние заряда. На основе данной 
модели проведен анализ движения потоков продуктов сгорания в области у передне- 
го днища и напряженно-деформированного состояния донного торца заряда. 


На функционирование РДТТ с зарядами из топлива с низкими значениями механиче- 
ских характеристик значительное влияние оказывают деформации зарядов под дейст- 
вием полетных перегрузок, перепадов давления и газодинамических нагрузок. В ряде 
случаев, например при высоких начальных температурах зарядов, использовании заря- 
дов большого относительного удлинения из низкомодульных топлив, нерациональных 
геометрических параметрах зарядов значительные локальные деформации зарядов и в 
первую очередь торцов зарядов могут являться причиной аномальной работы. Экспе- 
риментально установлено, что при нерациональном конструктивном оформлении РДТТ 
с высоким объемным заполнением в области переднего дна РДТТ, неэффективной схе- 
ме воспламенения существует вероятность перекрытия осевого канала заряда и отрыва 
донной манжеты от корпуса, что требует прогнозирования НДС зарядов в области пе- 
реднего торца. 

Определение взаимодействия между потоком продуктов сгорания и внутренней 
поверхностью конструкции требует расчета указанных процессов в нестационарной 
постановке. Для рассматриваемых объектов характерны такие особенности процессов, 
как значительные градиенты скорости и давления, наличие отрывных течений в соче- 
тании со сложностью геометрии полостей, наличие до-, транс- и сверхзвуковых зон в 
области течения. 

Для оценки ВБХ с учетом НДС зарядов разработана комплексная математическая 
модель, содержащая модели одномерного квазистационарного течения, двумерного те- 
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чения газов в области торцов зарядов, на основе которой рассчитано НДС зарядов в об- 
ластях, прилегающих к предсопловому и межсекционному объемам [1,2]. В работе на 
основе данной методики проведен анализ движения потоков продуктов сгорания в об- 
ласти у переднего дна и НДС донного торца заряда. 

Основу математической модели процессов нестационарного движения продуктов 
сгорания РДТТ составляют уравнения движения вязкого многокомпонентного газа. 
Система уравнений имеет вид [3]: 

- уравнение неразрывности 


Р ІЪ + рѴѴ % Га(І )> X е * = 7 -' 0; 

Ѵ1 к=і 

- уравнение количества движения 

с/ѴѴ „ 

Р~^ = -8гас1Р+ОіѵСТ-ХОі ѵ (т аі ІѴ ы ), Ѵх т еѴ о ,1>0; 


- уравнение энергии 


=-сІіѵ(РЧЧ)+сІіѵ(сг\и) + сііѵд, + М ■ ѴѴй* + 

+ X ^ іѵ ( СгѴѴо* ; - X рс к Е к ѴѴ 1>к ), \/х т &Ѵ а ,і> 0; 

к і 


(1) 


( 2 ) 


( 3 ) 


где Ѵ 0 - объем области; X - время, Е - ІД + \У 2 /2 - удельная полная энергия смеси, 

и = с ѵ Т ‘ удельная внутренняя энергия, ѴѴ - вектор скорости потока в данной точке 
(среднемассовая скорость); Р и Т - местные термодинамическое давление и температу- 
ра, х т - пространственные координаты; р - плотность среды; - вектор плотности те- 
плового потока; к=1,2,З...М - индекс компонента смеси; N - число компонентов смеси; 
Ск - массовые концентрации к-х компонентов, определяющие состав среды; сг - тензор 

напряжений вязкости; ѴѴ эк - вектор скорости диффузии компонента к (вектор скоро- 
сти компонента относительно потока: ѵѵ 0 к =ѵѵ к -ѵѵ, х - коэффициент теплопро- 
водности, с ѵ — удельная теплоемкость смеси, ц ѵ — удельное выделение теплоты, отне- 
сенное к единице массы, т вк - диффузионный поток массы к-ого компонента. 
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N 

- уравнение состояния: Р = рТ^ с Л. (4) 

к=1 

Для моделирования турбулентности используется стандартная двухпараметриче- 
ская модель, включающая уравнения для энергии турбулентности и скорости диссипа- 
ции. 

Для расчета рассматриваемых течений в каналах сложной формы, характерных 
для конструкций ЛА, использован модифицированный численный метод крупных час- 
тиц, реализованный на неравномерной сетке метода конечных элементов [4]. Введение 
указанной сетки объясняется необходимостью совместного моделирования процессов 
газодинамики и динамики деформирования заряда. 

Напряженно-деформированное состояние элементов РДТТ моделируется в дина- 
мической постановке с учетом нестационарного нагружения конструкции силовым 
воздействием со стороны продуктов сгорания ТТ. Для получения уравнения деформи- 
рования использован принцип виртуальных работ [4]: 

+ | 7 ]5и } с18 - ; ^ 

УѴ УѴ 5 уѵ 


где \ѵ и 8 объем и поверхность тела; сту, еу — компоненты тензоров напряжений и де- 
формаций, Ц| и а, - компоненты векторов перемещений и ускорений, и Т; — компонен- 
ты векторов объемных поверхностных сил, р - плотность материала. 

В области упругости связь напряжений и деформаций выражается в виде 

{ст} = [0 Е ].{е} , где [О ] - матрица упругих постоянных. В связи с тем, что рас- 
сматривается начальный период внутрибаллистического процесса, малый по продол- 
жительности, для определения сил сопротивления, вызванных явлениями вязкости, ис- 
пользовалась зависимость: 

= 2р(Т , 8 у , 8 у ){ё ц - д ] 8 п / З), ( 7 ) 


где Ц(Т, 8 у , 8у ) - коэффициент, учитывающий вязкость, 8у - символ Кронекера. 

Для описания процесса динамики деформирования приращение полной деформа- 
ции сіе у представляется в виде суммы [5] 

=Жг‘+<&; +<^ г , ( 8 ) 
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где бу - символ Кронекера, 8 е - приращение упругой деформации, е т - приращение де- 
формаций вследствие вязкости, г 1 - приращение температурной деформации. 

Для численного решения уравнения деформирования применялась итерационная 
процедура метода дополнительных деформаций. 

В соответствии с методом конечных элементов построена дискретная модель те- 
ла, состоящая из конечного числа связанных соответствующим образом в узловых точ- 
ках конечных элементов. Для расчета как механической, так и газодинамической задач 
в осесимметричной постановке использовались единые гибридные элементы, состоя- 
щие из нескольких кольцевых треугольных конечных элементов с тремя узлами с ли- 
нейной функцией распределения перемещений и температуры. 

Глобальное разрешающее уравнение МКЭ для НДС имеет вид 

[К К и) = !Р}- { В} т {|}-!Мг{^}, (9) 

где [К], {Р}, {О} , {М} — глобальные матрицы жесткости, узловых внешних сил, 
демпфирования и масс. 

Дискретизация по времени при численном расчете динамического процесса де- 
формирования осуществлена методом конечных разностей с применением неявной 
разностной схемы 


([К] + 


1 ЕО] + ^СЕМ]){и к »' ) = {Р к } + _1_{Г)> Т {ц 1 - 1 } + Л_( {М} Т ( 2 {и к }- {и ь -' }) . 


2М 


Дг 


( 10 ) 


По разработанным моделям создан алгоритм и программа численного моделиро- 
вания рассматриваемых процессов. Для оценки справедливости предложенных соот- 
ношений и алгоритма вычислений было произведено решения ряда тестовых задач и 
доказана сходимость и устойчивость решений. 

С использованием данной методики проведено решение задачи по оценке воздей- 
ствия потока продуктов сгорания заряда воспламенительного устройства, размещенно- 
го в переднем днище, на донный торец заряда, оптимизации параметров заряда воспла- 
менительного устройства и выбора рациональной конфигурации торца заряда. 

На рис. 1,2 представлена компьютерная визуализация полей давления и осевых 
напряжений в зоне донного торца заряда в начальный момент времени после срабаты- 
вания воспламенителя и графики радиальных перемещений точек, полученные при ис- 
следования взаимодействия продуктов сгорания воспламенителя РДТТ с торцом заря- 
да. 
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Давление 
6.7000е+06 

5.3800е+06 

4.0600е+06 

2. 7400е+06 
Щ 1.4200е+06 
™ 1.0000е+05 


Напряжения 5х 
Па 

3.3884е+06 

•2.1980е+08 

•7 7843е+06 

•1 3371 е+07 
Я -1 8957е-07 
™ -2.4543е+07 


Рис.1. Компьютерная визуализация полей давления и осевых напряжений в зоне 
переднего торца заряда при срабатывании воспламенителя (время процесса 
1=0.004 мс) 




Рис.2. Деформирование переднего торца заряда при срабатывании воспламените- 
ля (время процесса 1—0.004 с) и графики радиальных перемещений точек 1- 
3 торца 
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АННОТАЦИЯ. Изложен машинный эксперимент по изучению влияния механиче- 
ских характеристик твердого топлива и материалов корпусов маршевых ракетных 
двигателей на твердом топливе (РДТТ) при замене металлических корпусов на стек- 
лопластиковые и боралюмоуглепластиковые на боевую эффективность активно- 
реактивных снарядов (АРС) и надежность их функционирования. Даны статистиче- 
ский анализ работоспособности РДТТ и рекомендации по повышению технических 
характеристик изделий. 


Активно-реактивный снаряд (АРС) как техногенная система по принципу действия 
представляет собой систему, в которой сначала имеет место активное движение снаря- 
да и затем - активно-реактивное (рис. 1). 
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Рис. 1 . Декомпозиция системы АРС 
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Информационные данные о процессе отработки АРС калибра 1 52 и 202 мм с ме- 
таллическими корпусами свидетельствуют о наличии 3-5% отказов при проведении на- 
турных испытаний, проявляющихся в виде уменьшения дальности полета снаряда на 
10-12% или повреждений корпуса двигателя при движении на активном участке траек- 
тории. 

Анализ многоуровневой иерархической системы АРС (рис.1), показывает, что 
основная причина указанных отказов - потеря работоспособности его командного эле- 
мента - маршевого ракетного двигателя на твердом топливе (РДТТ). 

В процессе эксплуатации конструкция РДТТ испытывает динамические нагруз- 
ки как система с распределенными массами, обладающая реологическими деформаци- 
онными и прочностными характеристиками. При этом экстремальные нагрузки возни- 
кают при переходных режимах работы: стартовых перегрузках в период активного 
движения АРС, воспламенении твердого топлива (ТТ), выходе РДТТ на режим и т.п.. 
і.е. в периоды когда порядок эксплуатационных нагрузок определяется силами инер- 
ции. Колебательные процессы в двигателе, о чем свидетельствуют исследования внут- 

рикамерных (внутри корпуса РДТТ) явлений, изменяют характер физических процес- 
сов горения ТТ. 

Таким образом, можно предполагать, что потеря работоспособности двигателя 
связана с механической повреждаемостью, накапливаемой в его элементах при пере- 
ходных режимах эксплуатации и обусловлена динамическим характером нагружения. 

Механические повреждения в заряде ТТ приводят к аномалии его горения и, как 
следствие, к изменению вектора тяги, времени работы РДТТ и дальности полета АРС. 
а также к изменению характера физических процессов, приводящих к детонации топ- 
лива и разрушению корпуса АРС. 

Правомерность сделанных предположений подтверждается анализом прочности 
РДТТ при переходных режимах эксплуатации, а именно, на периоде активного движе- 
ния АРС, где имеют место указанные отказы[1].На рис. 2 представлены конструктив- 
ные схемы АРС калибра 1 52 и 202мм. 

Для указанных типов снарядов было выполнено исследование влияния механи- 
ческих характеристик ТТ и материалов корпусов АРС на прочность РДТТ [2], эксплуа- 
тационную надежность и боевую эффективность АРС. Оценка последних проводилась 
по дальности метания и дульной скорости, получаемых при использовании штатного 
"вышибного" заряда (без учета работы РДТТ). При этом определялась вероятность без- 
отказной работы РДТТ р(а‘"'\_ <[а]<'">,), где а<"> гаах <[<'>, - условие работоспо- 
собности двигателя; о ( "' ) , тах - максимальная интенсивность динамических напряжений 
в заряде ТТ; [о] ( '"Т - допускаемое напряжение при сжатии ТТ. 
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Рис. 2. Конструктивные схемы и варианты исполнения корпусов АРС калибра 
152мм (а ,б ,в) и 202мм (г) 

а ~ АРС со стальным корпусом: 1 -корпус, 2- сопловой блок, 3-ракетный 
двигатель, 4 -боевая часть, 5,6-прокладки;б -АРС с корпусом из компози- 
цинно-волокнистого материала (КВМ); в - АРС с составным по длине кор- 
пусом. 1 -металл , 2-КВМ; г — калибр 202 мм:1 -корпус, 2- сопловой блок, 3- 
ракетный двигатель,4-боевая 


В качестве материалов корпусов рассматривались сталь, стеклопластик и бора- 
люмоуглепластик (табл.1). Критерием выбора КВМ для изготовления корпуса служила 
длина снаряда: 4,6 калибра у 152мм АРС и 1 1,5 калибра у 202мм АРС. При этом для 
изготовления корпусов АРС калибра 152мм рассматривалось применение стеклопла- 
стика продольно-поперечного армирования (ППН 1:1) на основе волокна ВМ-1 и свя- 
зующего ЭХД, технологичного и сравнительно недорогого в изготовлении, имеющего 
достаточно высокие удельную прочность и жесткость. 
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Таблица 1. Механические характеристики материалов корпуса и двигателя АРС 


Материал 

Сталь 

Твердое топливо 

Стекло 

пластик 

ППН 1:1 

Боралюмо- 

углепла- 

стик 

ППН 1:2 



АРС 

152мм 

АРС 

202мм 

Модуль 
нормальной 
упругости Е , 

х ІО -2 МПа 

Ее 

2108 

6,25 

3,16 

450 

910 

Е г 

220 


Е х 

215 

420 

Коэффициент 
Пуассона р 

Ы 

0,31 

0,47 

0,49 

0,24 

0,16 

ВѲг 

0,21 

0,16 

йхѳ 

0,22 

0,14 


0,16 

018 

Ри 

0,14 

0,21 

Мит 

0,13 

0,21 

Допускаемое 

напряжение 

м?, 

хЮ' 2 МПа 

[ ст ѳ]„ 

5,7 

0,25 

0,19 

145 

9 

к]„ 

7,65 

6 

[ СТ Д, 

4,5 

4 

Объемный 

весу, кН/Ѵ 

У 

. 

78,2 

16,2 

17,5 

18,5 

21 


Здесь Е, ц - модуль нормальной упругости и коэффициент Пуассона, соответст- 
венно, изотропных материалов (корпуса, ТТ); Е Ѳ ,Е Г ,Е Х ^ 

[ а ѳ]р>[ ст г]р ’[ а лг]^~ модули нормальной упругости, коэффициенты Пуассона и преде- 
лы прочности при растяжении в тангенциальном Ѳ, радиальном г и осевом х направ- 
лениях, соответственно, материала корпуса АРС калибра 202мм (боралюмоуглепласти- 
ка). 

Для изготовления корпусов АРС калибра 202 мм, имеющих длину в 2,5 раза 
выше чем 152 мм АРС и вследствие этого обладающего пониженной изгибной жестко- 
стью, влияющей на точность поражения цели, рассматривалось применение боралюмо- 
углепластика (ППН 1 :2), в котором для продольного армирования используется бора- 
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люминиевый полуфабрикат, [3], а поперечное армирование осуществляется углерод- 
ным волокном ВМН-4. Такой материал практически не уступает стеклопластику по 
удельной прочности, но значительно превосходит его по удельной изгибной жесткости. 

Применение указанного композиционно-волокнистого материала при изготов- 
лении несущего кожуха опытного ствола 30мм артиллерийского автомата позволило 
повысить технические характеристики ствола, влияющие на точность стрельбы, что 
подтверждено результатами стендовых и полигонных испытаний [4]. 

Из рис.З видно, что кривая ДА - Ад имеет экстремум, что позволяет сделать 
предположение о существовании оптимальной массы (^ = 32ка), обеспечивающей мак- 
симальное увеличение дшіьности метания АРС калибра 152мм на 12,1% (опытный 2- 
рис.2в, таол.2) вследствие частичной замены стального корпуса штатного снаряда 
стеклопластиковым (табл.1, рис.2в). В то же время указанный вариант исполнения 
обеспечивает значительное увеличение АѴ д ,% (на 24,9%), что значительно увеличивает 

могущество боеприпаса и его боевую эффективность артиллерийского и ракетного 
оружия. 

Полная замена металла в корпусе на стеклопластик (опытный 1 - рис.2б) приво- 
дит к значительному снижению массы Д^,% = 58,1% и увеличению дульной скоро- 
сти АѴ д ,% =56,1% снаряда. Однако, имеющее место снижение инерционности АРС, об- 
ладающих малой массой, вследствие расходования кинетической энергии на преодо- 
ление сил сопротивления, сокращает дальность метания (ДА иеш ,%= 26,4%). В то же 
время, применение АРС с указанным корпусом может оказаться эффективным при 
стрельбе на малые расстояния, т.к. высокая дульная (начальная) скорость способствует 
высокой пробиваемости цели. 



Рис .3. Зависимость изменения дальности метания Аз мст от изменения Аа массы 
АРС 
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Рис. 4. Зависимость дульной скорости Ѵ д от массы с/ 152мм АРС 

В то же время полная или частичная замена в корпусе снаряда металла на стек- 
лопластик приводит к некоторому повышению максимального уровня интенсивности 
напряжений в ТТ а,^ по причине увеличения осевых а х и центробежных ха х пере- 
грузок (табл. 2), что приводит к незначительному снижению вероятности безотказной 
работы РДТТ. Повышение вероятности безотказной работы двигателя АРС калибра 
152мм может быть достигнуто путем внесения конструктивных изменений в РДТТ, на- 
пример изменением количества и толщины прокладок поз. 5,6 (рис.2а), приводящее к 
снижению максимального уровня напряжений в ТТ [2]. Замена стального корпуса 
202мм АРС на боралюмоуглепластиковый приводит к снижению массы АРС на 35,9%, 
повышению дульной скорости на 27,2% и дальности метания на 29,2% (табл.2). Имею- 
щее при этом место, снижение максимального уровня интенсивности напряжений в ТТ 
связано с увеличением перегрузок, приводящих к повышению контактных давлений на 
границе корпус РДТТ -ТТ (опытный 1 и опытный 2, табл. 2). Это, в свою очередь, по- 
вышает вероятность безотказной работы двигателя и эксплуатационную надежность 
АРС. Увеличение дальности метания снаряда калибра 202мм позволяет, уменьшив мас- 
су заряда ТТ, величина объемного веса которого практически не отличается от объем- 
ною веса боевого заряда (ВВ), за счет его длины, например в два раза ( Д І іг = -50 % . 
табл. 2 -опытный 2), повысить могущество АРС, увеличив массу ВВ кратюГвышеука- 
занному. Уменьшение длины заряда ТТ за счет улучшения контакта с корпусом приве- 
дет к еще большему по сравнению со штатным снарядом снижению максимального 

уровня интенсивности напряжений в ТТ (а^\ тах /[о]Ч=0,51), повышению вероятно- 
сти без- отказной работы (до 0,91) и, как следствие, повышению эксплуатационной на- 
дежности РДТТ и всего АРС в целом. 
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Таблица 2. Конструктивные, динамические, прочностные характеристики и пара- 
метры движения АРС 


№ 

п/п 

Характеристики 
и параметры 

АРС 

Калибр АРС 

152мм 

202мм 

штат- 

ный 

опыт- 
ный 1 

опыт- 
ный 2 

штат- 

ный 

опыт- 
ный 1 

опыт- 
ный 2 

1 

Я, кг 

42,7 

18,2 

32 

135,7 

87 

87 

2 

^ СН ’ ММ 

700 

700 

700 

2332 

2332 

2332 

3 

Ь дв ,мм 

170 

170 

170 

1356 

1356 

678 

4 

а х ,§ • 10 ~ 3 м/с 2 

10,74 

24,42 

19,54 

1,26 

2,47 

2,47 

5 

та х , рад / с 

1062 

1579 

1210 

53,2 

61 

61 

6 

1 д ,мс 

14,4 

9,6 

10,5 

70,3 

36,3 

36,3 

7 

Ѵ д ,м/ с 

770 

1202 

962 

297 

378 

378 

8 

5 м Л0~\м 

14,053 

10,336 

15,753 

6,592 

8,520 

8,520 

9 

° {т) . т ,МПа 

24,3 

24,6 

24,45 

18,6 

17,8 

9,8 

10 

[о] {т) с ,МПа 

25 

25 

25 

19 

19 

19 

11 

АЯ,% 

0 

-58,1 

-25,6 

0 

-35,9 

-35,9 

12 

АѴ д ,% 

0 

56,1 

24,9 

0 

27,2 

27,2 

13 

А Ь дв ,% 

— 

— 

— 

— 

— 

-50 

14 

Ь$мет ’ 0// ° 

0 

-26,4 

12,1 

0 

29,2 

29,2 

15 

Да,- ,% 

чпах 

0 

1,23 

0,62 

0 

-4,3 

-47,3 

16 

«Н- /мН 

0,972 

0,984 

0,978 

0,98 

0,94 

0,51 

17 

Р(° (т \„ <[сг] (т) с) 

0,70 

0,72 

0,71 

0,86 

0,88 

0,91 


ВЫВОДЫ 

1. Полная или частичная замена металлического корпуса активно-реактивного 
снаряда на композиционно волокнистый может способствовать повышению его балли- 
стических характеристик: увеличению дальности поражения и скорости встречи снаря- 
да с целью и, как следствие - повышению боевой эффективности и могущества ракет- 
но-артиллерийского оружия. 
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2. Использование боралюминиевого полуфабриката для продольного армирова- 
ния в композиции с углеродным волокном повышает изгибную жесткость корпуса сна- 
ряда, влияющую на точность поражения цели, работоспособность РДТТ и эксплуата- 
ционную надежность АРС. 
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АННОТАЦИЯ. Разработана методика компьютерного моделирования взаимодейст- 
вия ударников с преградами из различных конструкционных материалов. Спрогно- 
зированы результаты взаимодействия, а также выявлено влияние материалов сер- 
дечников на пробивное действие ударников. Установлено, что за счет изменения 
компоновки и формы ударника для дозвуковой скорости взаимодействия можно су- 
щественно увеличить его пробивное действие. 


В настоящей работе при помощи методов математического моделирования [1,2] иссле- 
дуется процесс пробития преград осесимметричными комбинированными ударниками 
по нормали. Рассматриваются ударники оживальной и затупленной формы, оболочеч- 
ные и полуоболочечные. Конструктивно ударники являются телами вращения и состо- 
ят из трех элементов: стальной оболочки, свинцовой рубашкой и сердечника. Диаметры 
всех составных частей ударников равны. Внешний диаметр оболочки равен 9,25 мм, 
толщина - 0,5 мм, масса - 3,53 г, материал - мягкая сталь. Диаметр сердечника равен 
7,4 мм. Масса рубашки - 1,6 г. Материал сердечников: сталь У10А, сталь 10, свинец, 
вольфрамовый сплав (ВНЖ), уран, золото. Материал преград: сплав Д16, сталь 3, высо- 
копрочная сталь (ВПС), титан. Скорость удара - 290 м/с. 

Экспериментальные исследования проведены на баллистическом стенде с ис- 
пользованием пороховой метательной установки [3]. Начальные и запреградные скоро- 
сти ударников регистрировались с помощью индукционного измерителя и рам- 
мишеней. Ошибка измерителя скорости не превышала 1%. В экспериментах использо- 
ваны ударники с сердечниками из стали У 10 А, преграды из сплава Д16 толщиной 4 мм 
и 8 мм и преграды из стали 3 толщиной 4 мм. На рис. Іа приведены обозначения и схе- 
мы ударников в разрезе. Здесь 1- сердечник, 2 - оболочка, 3 - рубашка. Ударники ожи- 
вальной формы обозначены как В1, В2, ВЗ. Затупленной формы - В4. Ударник В2 - 
полуоболочечный, а ударник ВЗ — оживальный с затупленным сердечником. Массы со- 
ставных элементов всех ударников, использованных в эксперименте равны. Результаты 
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измерений запреградной скорости, пересчитанные для начальной скорости удара 290 
м/с приведены в таблице 1 

Детальный анализ процесса пробития преград проведен при помощи методики 
компьютерного моделирования, базирующейся на основных положениях механики 
сплошных сред в рамках двумерной постановки для осевой симметрии в лагранжевых 
независимых переменных [1,2]. Разработанный программный комплекс позволяет про- 
гнозировать действие ударного и взрывного нагружения материала с учетом волновых 
эффектов и фрагментарного разрушения при реализации отрывного или сдвигового ме- 
ханизмов, что дает возможность избежать нефизических деформаций и адекватно опи- 
сать таким образом изучаемые процессы и явления. Необходимые физико- 
механические характеристики рассматриваемых материалов взяты из работ [4,5]. 
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Результаты вычислений в виде текущих конфигураций ударник-мишень пред- 
ставлены на рис. 16 для преград из сплава Д16 толщиной 8 мм, на рис.ів для сплава Д16 
толщиной 4 мм и на рис.іг для стали 3 толщиной 4 мм. Видно, что во всех случаях, как 
и в эксперименте [3] пробитие преград для ударников В1 и В2 происходит по механиз- 
му «прокола», а для ударников ВЗ и В4 по механизму срезания «пробки». Форма и ос- 
таточные размеры оболочки, сердечника и пробоины, полученные расчетным и экспе- 
риментальным путем, соответствуют друг другу. На рис. 16 видно, что для данного ва- 
рианта имеет место отскок оболочки от лицевой поверхности преграды, обусловленный 
упругой реакцией материала. 

В таблице 1 приведены рассчитанные значения запреградных скоростей центров 
масс ударников. Хорошее совпадение с экспериментом получено в случае стальной 
преграды, расхождение не превышает 2,7%. Для преград из сплава Д16 наибольшее 
расхождение расчетных и экспериментальных данных достигает 4,6% (толщина 4 мм, 
ударник ВЗ) и 3,4% (толщина 8 мм, ударник ВЗ). 

На рис.2а представлены рассчитанные зависимости скорости центра масс сер- 
дечников от времени при пробитии преграды из сплава Д16 толщиной 4 мм. Вид этих 
зависимостей отражает характерные для каждого типа ударников способы пробития 
преграды. Так для ударников В1 и В2 кривые скорости центра масс вначале более по- 
логи, чем аналогичные участки кривых для типов ВЗ и В4. Это вызвано тем, что сила 
сопротивления внедрению ударников с оживальной головной частью возрастает плав- 
но вслед за плавным увеличением поперечного сечения сердечника в отличие от удар- 
ников с затупленной головной частью. При внедрении ударников с оживальной голов- 
ной частью деформируется и уплотняется объем материала преграды вокруг оживала 
по всей его длине, а при внедрении ударников с затупленной головной частью проис- 
ходит сдвиг материала преграды, приводящий к образованию пробки по диаметру сер- 
дечника. к 

Таблица 1. 


Преграда 

Д16 
(4 мм) 

Д16 

(8мм) 

Ст.З 
(4 мм) 

Ударник 


В 1 
(м/с) 

Расчет 

196 

98 

167 

Эксп - т 

190 

95 

166 

В2 

(м/с) 

Расчет 

232 

131 

182 

Эксп - т 

241 

135 

186 

ВЗ 

(м/с) 

Расчет 

226 

90 

149 

Эксп - т 

216 

87 

147 

В4 

(м/с) 

Расчет 

253 

174 

183 

Эксп - т 

- . 

258 

173 

188 
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Рис. 2. 

Далее с целью прогноза смоделирован процесс взаимодействия показанных на 
рис. Іа ударников с 4-х мм преградами из ВПС (динамический предел текучести равен 
0,9 Гпа) и титана. Рассчитанные конфигурации ударник-мишень для преград из ВПС 
приведены на рис. За. Видно, что в этом случае для ударников В2 и В4 имеет место от- 
скок оболочки от преграды, а для ударников В1 и ВЗ оболочки деформируясь сминают- 
ся на лицевой поверхности преграды. Ударник В1 пробивает преграду и движется 
дальше со скоростью 52 м/с. Снижение скорости ударника В2 заканчивается через 160 

мкс от начала процесса, ударника ВЗ - через 180 мкс, ударника В4 - через 60 мкс 
(рис.2б). 

Рассчитанные конфигурации ударник-мишень для варианта с титановой прегра- 
дой приведены на рис.Зб. Видно, что динамика пробития и картины разрушения пре- 
град и ударников мало отличаются от вариантов расчетов, представленных на рис.ів и 
рис. 1г. Снижение скорости ударника В1 в этом случае заканчивается через 270 мкс от 

начала процесса, ударника В2 - через 1 10 мкс, ударника ВЗ - через 200 мкс и ударника 
В4 - через 60 мкс. 

Вычисленные значения запреградных скоростей центров масс ударников приве- 
дены в табл.2. Видно, что также как и для преград из сплава Д16 наибольшим пробив- 
ным действием обладает ударник В4, а наименьшим - В1. Причем в вариантах с тита- 
новой преградой для ударников В1 и В2 запреградные скорости отличаются на 46%, а 
для ударников ВЗ и В4 - на 5,4%. 
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Рис.З. 


Таблица 2. 


Уд - ник 

Д16 

(м/с) 

Ст.З 

(м/с) 

впс 

(м/с) 

Титан 

(м/с) 

В1 

196 

167 

52 

89 

В2 

232 

182 

83 

130 

ВЗ 

226 

149 

129 

185 

В4 

253 

183 

151 

195 


Для выявления влияния материала сердечника на пробивное действие ударника 
смоделирован процесс взаимодействия ударников В 1 и В4 с 4-х мм преградой из стали 
3. В качестве материалов сердечников рассмотрены сталь 10, свинец, уран, ВНЖ и зо- 
лото. Выбор последних трех материалов продиктован соображением целесообразности 
сравнения механизмов действия ударников с тяжелыми прочными (ВНЖ) и слабопроч- 
ными (Аи) сердечниками. Рассчитанные текущие конфигурации ударник-мишень для 
ударников с сердечниками из ВНЖ, урана и золота представлены на рис. 4. 

В табл.З приведены рассчитанные значения запреградных скоростей ударников с 
сердечниками из шести различных материалов. Во второй строке таблицы приведены 
массы ударников. Получено, что запреградная скорость ударника с сердечником из 
стали 10 в несколько раз меньше, чем для аналогичного ударника с сердечником из 
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В1) 
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Рис.4. 
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Таблица 3. 


Сердеч- 

к 

У10А 

Ст.10 

РЬ 

ВНЖ 

И 

Аи 

Уд-ник 

Масса 

(г) 

14,96 

14,96 

19,3 

26,73 

28,46 

29,33 

В1(м/с) 

167 

24 

0 

233 

241 

163 

В4 (м/с) 

185 

50 

0 

249 

253 

197 


У10А. Ударник со свинцовым сердечником растекается в радиальном направлении по 
преграде, образуя в ней вмятину, но, не пробивая ее. Значение запреградных скоростей 
для ударников с золотыми сердечниками получились близкими к значениям для удар- 
ников с сердечниками из стали У10А, несмотря на то, что к концу процесса золотые 
сердечники были сильно деформированными и имели грушевидную форму (рис.4в). 
Ударники с сердечниками из ВНЖ (рис. 4а) имели запреградную скорость на 39% и 
32% выше, чем ударники с сердечниками из стали У10А, хотя масса последних в 1,8 
раза меньше Самое высокое значение запреградной скорости равное 253 м/с получено 
для ударников с сердечниками из урана (рис. 46). В этом случае для ударника В1 на- 
блюдается отскок оболочки от преграды, пробитие проходило по механизму “прокола”, 
снижение скорости закончилось к 180 мкс. Для ударника В4 пробитие проходило по 
механизму срезания “пробки”, а снижение скорости закончилось к 40-й мкс. 

Таким образом, установлено, что за счет изменения компоновки и формы ударника 
для дозвуковой скорости взаимодействия можно существенно увеличить его пробивное 
действие. Кроме этого, при помощи разработанной методики компьютерного модели- 
рования спрогнозированы результаты взаимодействия рассмотренных ударников с пре- 
градами из различных конструкционных материалов, а также выявлено влияние мате- 
риалов сердечников на пробивное действие ударников. 
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АННОТАЦИЯ. Приведены результаты экспериментальных исследований поведе- 
ний углерод-углеродных материалов при ударно-волновом нагружении. Получен- 
ные результаты показывают высокую чувствительность структуры фронта ударных 
волн и откольной прочности к структуре композиционных материалов. 


Композитные материалы на основе углеродных волокон находят широкое применение 
в качестве элементов конструкций, испытывающих высокие силовые и тепловые на- 
грузки. Поэтому прогнозирование результатов импульсного воздействия на эти конст- 
рукции является важной практической задачей, для решения которой требуется знание 
реологических и теплофизических свойств вещества. Однако имеющийся набор опыт- 
ных данных о поведении композитных материалов при ударно-волновом воздействии 
крайне малочисленен [1-4]. В данной работе проведены эксперименты по исследова- 
нию структуры фронта ударных волн и откольной прочности в трех различных угле- 
род-углеродных материалах (УУМ) при давлении порядка 2 ГПа (20 кбар). 

^ УМ представляют собой объемно-армированные композиты, отличающиеся 
между собой способом плетения углеродных волокон в пирографитовой матрице. 

Схема экспериментов приведена на рис. 2. С помощью плосковолнового взрывно- 
го генератора 1 в стальном экране 2 диаметром 120 мм и толщиной 20 мм формирова- 
лась ударная волна амплитудой 9 ГПа и длительностью около 10 мкс. Исследуемые об- 
разцы 3 имели размер 60x60 мм, толщину 5-10 мм, находились в контакте с экраном и 
были ориентированы таким образом, чтобы плоскость расположения армирующих уг- 
леродных волокон была параллельна фронту ударной волны. Структура фронта удар- 
ной волны и величина откольной прочности определялись по скорости свободной по- 
верхности. Для регистрации скорости использовался лазерный интерферометр ѴІ8АК 
[5] с постоянной интерферометра, равной 80,8 м/с, что позволило провести измерения с 
точностью ± 2 м/с и временным разрешением около 3 нс. Для отражения лазерного лу 
ча, на поверхность образца наклеивалась алюминиевая фольга 4 толщиной 7 мкм. Луч 
интерферометра фокусировался в пятно диаметром около 100 мкм. 
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Рис. 1. Структура образцов УУМ 
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Рис. 2. Схема экспериментов и профиль скорости свободной поверхности экрана 
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Для определения параметров ударной волны, входящей в образец, предваритель- 
но проведены измерения скорости свободной поверхности стального экрана при отсут- 
ствии образца (рис.2). Видно, что в экране формируется двухволновая структура, обу- 
словленная упруго-пластическими свойствами стали (СтЗ), поэтому сначала на свобод- 
ную поверхость выходит упругий предвестник (амплитуда которого равна 1.3 ГПа), а 
примерно через 0.6 мкс пластическая волна сжатия. Если между экраном и образцом 
обеспечить акустический контакт, то такая двухволновая конфигурация будет входить 
в образец. В действительности же, из-за наличия воздушной прослойки, упругий пред- 
вестник затухает и в начальный момент времени в образце формируется ударный ска- 
чок с характерной шириной около 100 нс. На рис. 2 видно также, что после достижения 
максимума, скорость свободной поверхности экрана начинает уменьшаться, что обу- 
словлено волной разгрузки со стороны плосковолнового генератора. В результате, в 
экране возникают растягивающие напряжения, что приводит к разрушению - отколу, 
проявляющемуся на профиле скорости свободной поверхности в виде характерного от- 
кольного импульса и последующих колебаний скорости. Величина откольной прочно- 
сти СтЗ, найденная из данного опыта, равна 1.5 ГПа. 

Результаты экспериментов с УУМ приведены на рис.З. 
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Там же указаны начальная плотность материала р 0 и скорость звука с 0 , измерен- 
ная на ультразвуковой установке с использованием сигналов высокой амплитуды (~1 
бар). Для мелко- и среднеструктурного материалов скорость звука имеет величину, ти- 
пичную для композитов подобного типа, около 3 км/с, тогда как в крупноструктурном 
материале зарегистрирована аномально высокая скорость звука - почти 7 км/с, что го- 
раздо больше, чем у всех подобных материалов. Одна из гипотез, объясняющих столь 
высокую скорость распространения возмущений - это то, что малые возмущения рас- 
пространяются вдоль углеволокна, имеющего, как известно, очень высокий модуль уп- 
ругости. Однако даже в этом случае дополнительно необходимо предположить, что 
значительная доля волокон ориентирована под большими углами к плоскости образца. 

Отличительной чертой всех трех профилей скорости является наличие предвест- 
ника, амплитуда которого монотонно возрастает от нуля до максимального значения 
равного примерно 0.2 ГПа для мелко- и среденеструктурного, и 0.5 ГПа для крупно- 
структурного. Причем, его наличие, как уже отмечалось, не связано с упругим пред- 
вестником в стальном экране, а обусловлено механическими свойствами композицион- 
ных материалов. Вслед за предвестником регистрируется относительно резкое нараста- 
ние скорости до максимального значения в мелко- и среднеструктурном материалах. В 
материале с крупной структурой не удалось записать полностью фронт ударной волны, 
поскольку возмущения плоскостности волны за счет структуры материала или его раз- 
рушения на фронте оказались слишком высоки. Характерная ширина фронта ударной 
волны при заданных условиях нагружения определяется структурой вещества и, как 
видно из рис.З, мелкая и средняя структуры в этом смысле существенно отличаются 
друг от друга: ширина фронта для средней структуры превышает 1 мкс, тогда как для 
мелкой она примерно вдвое меньше. Столь заметное различие является следствием то- 
го, что мелкоструктурный материал имеет отчетливо выраженную продольно- 
ориентированную структуру, тогда как среднеструктурный - более поперечно- 
ориентирован. Если для описания эволюции волн малой амплитуды в этих материалах 
применять модель вязкой жидкости, то их коэфициенты вязкости также будут отли- 
чаться вдвое [6]. Вероятно, более корректными для описания подобных материалов яв- 
ляются релаксационные модели [7], основной параметр которых - время релаксации - 
также извлекается из анализа структуры фронта ударной волны. 

Выполненные оценки являются корректными лишь при выполнении условия ста- 
ционарности ударных волн, что вообще говоря требует дополнительной проверки, в 
частности, проведения измерений при изменении толщины образцов. Это обстоятель- 
ство не является принципиальным для оценки ширины фронта, но без подобных экспе- 
риментов сложно обсуждать природу предвестника, отчетливо регистрируемого для 
всех трех композитов. Если, например, он «убегает» вперед, то следует говорить о на- 
личии особенности на ударной адиабате, аналогичной пределу текучести для упруго- 
пластических сред. Причем ответить на вопрос, будет ли этот предел действительно 
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связан с обратимым переходом материала в область пластического течения или он сов- 
падает с порогом разрушения на сжатие, можно лишь в результате выполнения доста- 
точно сложных и многочисленных экспериментов. Если же предвестник образует еди- 
ный стационарный комплекс, то его формирование обусловлено особенностях кинети- 
ки релаксационных процессов во фронте ударной волны, которые также могут носить 
как обратимый, так и необратимый характер. 

Проведенные эксперименты наряду с шириной фронта ударной волны позволяют 
определить величину откольной прочности в средне- и мелкоструктурном материале. 
Поскольку входящие в образцы импульсы сжатия являлись треугольными (рис.1), то их 
взамодействие со свободной поверхностью приводит к растягивающим напряжениям, 
максимальное значение которых (откольная прочность) определялось по формуле: 
<7=0.5р 0 СоД\Ѵ - где Д\Ѵ - перепад скорости свободной поверхности между максимумом 
и ее значением перед откольным импульсом (рис.З). Любопытно отметить, что, как и 
ширина фронта ударной волны, откольная прочность среднеструктурного материала в 
два раза больше, чем мелкоструктурного. Это может быть отражением того факта, что 
скорость образования несплошностей при импульсном растяжении лимитируется ско- 
ростью вязкого роста пор [8] и если предположить, что ширина фронта ударной волны 
и скорость роста пор определяются одной и той же вязкостью, то полученный резуль- 
тат является очевидным. Отметим, что материалы сохранили, по крайней мере частич- 
но, прочность на сжатие. 

Для дополнительного анализа механического поведения материалов также были 
проведены измерения диаграмм деформирования при статическом сжатии. Результаты 
представлены на рис. 4. 

Отметим, что поскольку сжатие происходило параллельно плоскости армирова- 
ния, то заметной поперечной деформации не наблюдалось (коэффициент Пуассона 
близок к нулю). Сжатие проводилось до появления первых признаков разрушения. 
Крупноструктурный образец разрушился до отдельных прядей и рассыпался, на сред- 
неструктурном было заметно расслоение между слоями армирования, а мелкострук- 
турный сохранился практически без изменений, что наблюдается и по кривым разгруз- 
ки. Материалы оказались сравнительно «мягкими» по сравнению с монолитными мате- 
риалами - модуль Юнга составил порядка 0,7 ГПа и скорость звука, вычисленная по 

соотношению с = у[Ё[р получилась на уровне 0,7 км/с, что типично для высокопорис- 
тых материалов. Отклонения на начальных участках и в начале кривых разгрузки свя- 
заны с особенностями испытательной машины; поведение материалов при сжатии 
практически линейно вплоть до нагрузок порядка 1-1,5 кбар и очень высокой деформа- 
ции - 20-25%, причём мелкоструктурный материал практически не потерял сопротив- 
ления на сжатие. 
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Рис. 4. Диаграммы статического сжатия 

Сопоставление ударноволновых и статических данных по сжатию материалов го- 
ворит о том, что наблюдаемый на фронте ударной волны предвестник действительно 
имеет упругий характер, поскольку лежит в том же диапазоне амплитуд, что и статиче- 
ские пороги начала разрушения - по-видимому, происходит разрушение пироуглерод- 
ной матрицы; однако в среднеструктурном материале, а в особенности мелкоструктур- 
ном, поперечно ориентированные волокна остаются неразрушенными и при амплиту- 
дах волны ~ 20 кбар. Крупная же структура плетения волокна не препятствует разру- 
шению, причем за счёт связующего, а не волокна, которое в условиях сжатия, как ста- 
тического, так и ударно-волнового, будет подвергаться большему деформированию. 

Таким образом, проведенные эксперименты показывают высокую чувствитель- 
ность структуры фронта ударных волн и откольной прочности к структуре композит- 
ных материалов, а подходы механики сплошных сред, разработанные для описания мо- 
нолитных материалов при ударно-волновых воздействиях, применительно к композит- 
ным, а в особенности к пористым, нуждаются в серьёзной доработке. 
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АННОТАЦИЯ, Приводится математическое описание модели термоэррозионной и 
механической стойкости стволов артиллерийских автоматов в условиях стрельбы 
короткими или длинными очередями с темпом стрельбы 600 выстр./мин. Комплекс- 
ная математическая модель включает описание процессов внутренней баллистики, 
теплообмена, теплопроводности и износа с учетом структурно-фазовых превраще- 
ний в материале орудийной стали, а также механических и температурных напря- 
жений 1-го и 2-го рода, возникающих при циклическом высокоимпульсном тепло- 
вом нагружении. Результаты численных исследований по разработанной модели 
иллюстрируются примером расчетов 30-мм гипотетической автоматической пушки, 
имеющей начальную скорость снаряда 960 м/с и максимальное давление 380 МПа. 
Ствол изготовлен из орудийной стали 38ХНМФА. 


Существующий в проектной практике метод расчета стволов стрелкового и артилле- 
рийского оружия на прочность основан на макроструктурном подходе, основанном на 
гипотезе о том, что живучесть ствола определяется снижением только его баллистиче- 
ских характеристик вследствие износа поверхности канала ствола под действием теп- 
ловых, силовых нагрузок, а также химического взаимодействия поверхностного слоя 
металла с продуктами сгорания пороховых газов. При этом интенсивность износа зави- 
сит от локальных характеристик газового потока и тепловых и механических характе- 
ристик металла. 1 акой традиционный подход не учитывает микроструктуру металла по 
его толщине и ее изменение в результате циклического теплового и силового воздейст- 
вия и ее изменение в результате структурно-фазовых превращений и, как показывает 
опыт, дает завышенные значения циклической долговечности ствола [1]. 

При макроструктурном подходе к оценке прочности материал стенки ствола 
представляется как сплошная среда со среднеинтегральными характеристиками проч- 
ности и пластичности, определяемыми при его механических испытаниях в статиче- 
ских условиях. Оценка работоспособности конструкции выполняется на основе анализа 
ее напряженно-деформированного состояния, компоненты которого получаются на ос- 
нове указанного подхода. Так, в стволе - толстостенной трубе, механические напряже- 
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ния, обусловленные действием давления пороховых газов, достигают экстремальных 
значений (Мизес - Р, рис Л) на его внутренней поверхности, что и определяет положе- 
ние опасных точек. 

Накопление повреждений в результате знакопеременных циклических напряже- 
ний приводит к снижению долговечности и механической живучести. Следует отме- 
тить, что кроме рассматриваемых термических напряжений, связанных с температур- 
ными градиентами, и называемых напряжениями первого рода, в металлах, вследствие 
их неоднородности, обусловленной поликристаллическим строением, возникают тер- 
мические микронапряжения, называемые напряжениями второго рода или термострук- 
турными. Результаты экспериментального исследования стали 12Х18Н10Т на образцах 
при изотермическом термоциклировании, т. е. при отсутствии термонапряжений перво- 
го рода, в интервале температур 300 - 900 К свидетельствуют о присутствии в металле 
термоструктурных напряжений высокого уровня, подтверждением чему служит боль- 
шое количество повреждений на поверхности образцов: микротрещин - с 1500 термо- 
циклов (тц), и макротрещин - с 5000 тц (рис. 2) [2]. 

Природа термических напряжений, связанных с теплосменами в результате че- 
редования термического воздействия газового потока и охлаждения стенки, обусловле- 
на градиентами температуры по её толщине. При этом в наиболее нагретых слоях 
стенки, расположенных в непосредственной близости к поверхности канала ствола, 
действуют сжимающие, а в менее нагретых (внутри стенки) - растягивающие напряже- 
ния, т. е. имеет место разгрузка поверхности канала ствола со смещением опасных то- 
чек в радиальном направлении, что способствует повышению его несущей способности 
(Мизес -РТ, рис. 1). 


700 



X МизесК 


Мизес-Р 

МизесЖ 


Мизес-РТ 


Ж : 1 і- 

1.5 2,0 2,5 3,0 3,5 4,0 

Радиус,см 


4,5 


Рис. 1 . Интенсивности напряжений по критерию Мизеса для гладкостенного 
ствола 


ХИМИЧЕСКАЯ ФИЗИКА И МЕЗОСКОПИЯ. Том 7, № 3 


267 


Д.П. АЛЕКСАНДРОВ, В.Ф. ЗАХАРЕНКОВ, В С. ЦЕПЕЛЕВ 



Рис. 2. К экспериментальному исследованию термоструктурных напряжений при 
изотермическом нагреве 

а - схема установки для термоусталостных испытаний при нестесненных 
температурных деформациях: 

1 - основание; 2 - ванна холодильника; 3 - испытываемый образец; 4 - 
электропечь; 5 - рама; 6 - роликовая направляющая рамы; 7 - подъемный 
трос; 8 - электродвигатель подъемно-транспортного механизма рамы с об- 
разцом; 9 - блок управления; 
б - образец для испытаний; 

в - поверхность образца после 2500 теплосмен (увеличение в 100 раз) 

Расчетные макротермические напряжения, значения которых обычно исполь- 
зуются при оценке прочности стволов артиллерийских орудий, снижают уровень ин- 
тенсивности напряжений в опасных точках на поверхности канала ствола во время вы- 
стрела и поэтому приводят к заниженной оценке механического износа ствола. Однако 
в реальных условиях в стенке нагретого ствола одновременно с указанными напряже- 
ниями действуют и так называемые термоструктурные напряжения, обусловленные 
различием коэффициентов термической деформации (коэффициентов линейного рас- 
ширения) главных структурных компонентов легированной стали (железа и карбидов 
металлов), а также образованием взаимно переходящих друг в друга структур мартен- 
ситного и аустенитного превращения. 

Расчет температурных полей по толщине стенки ствола и износа проводился с 
помощью программного комплекса СТРИЖ, включащего программы расчета внутрен- 
ней баллистики [3], теплообмена [4], теплопроводности, диффузии и термоизноса [5]. 
На рис. 3 показан характер изменения мартенситного и аустенитного превращения в 
различных слоях по толщине стенки ствола в поперечном сечении, соответствующем 
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Рис. 3. Характер изменения мартенситного и аустенитного превращения в ору- 
дийной стали 38ХНМФА при непрерывной стрельбе очередью 120 выстре- 
лов с темпом 10 выс/с: 0 — изменение относительного содержания мартен- 
сит на поверхности канала ствола; 1 — изменение относительного содер- 
жания мартенсита в слое на глубине 50 мк; 2 — изменение относительного 
содержания мартенсита в слое на глубине 75 мк; 3 — изменение относи- 
тельного содержания мартенсита в слое на глубине 1 00 мк; 4 — изменение 
относительного содержания мартенсита в слое на глубине 125 мк; 5 - из- 
менение относительного содержания мартенсита в слое на глубине 150 мк; 

6 - аустенитное превращение на поверхности канала ствола с учетом изно- 
са 


приложению максимальной газодинамической и термомеханической от выстрела (5 ка- 
либров от начала калиберной части). 

Рассмотренные режимы стрельбы приводят к термоэрозионнму износу поверх- 
ности канала ствола. Характер изменения радиального износа и максимальной скоро- 
сти износа в зависимости от текущего выстрела для непрерывной очереди представлен 
на рис. 4 и рис. 5. 

Полученные данные позволяют перейти к оценке механической живучести ство- 
лов при импульсном высокотемпературном нагружении. 
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Рис.4. Характер изменения термического износа от числа выстрелов в очереди 



Номер выстрела 


Рис.5. Характер изменения максимальной скорости термического износа в зави- 
симости от номера выстрела в непрерывной очереди 

Наличие в слоях ствола различных линейных температурных деформаций при- 
водит к появлению дополнительных напряжений растяжения-сжатия и, как следствие 
этого, к снижению циклической долговечности ствола и уменьшению его механиче- 
ской живучести. Признаком таких напряжений при циклическом термическом и сило- 
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вом нагружении является постепенное появление и развитие сетки трещин на внутрен- 
ней поверхности канала ствола. Для решения проблемы формирования таких трещин 
необходимо знать уровень циклических напряжений в режиме нагрев-охлаждение. Как 
показали расчеты, термоструктурные напряжения между выстрелами и к началу теку- 
щего выстрела могут достигать предела текучести материала ствола, а непосредственно 
при выстреле даже превосходить его. 

Для оценки уровня термических напряжений в легированной стали при изотер- 
мическом нагреве вследствие различия коэффициентов термического расширения меж- 
ду фазовыми составляющими предлагается структурная модель указанной стали в виде 
дисперсно-упрочненного композита (КМ) с матрицей из углеродистой стали (ж) и кар- 
бидами легирующих элементов (к) в качестве наполнителя (рис. 6) [6], [7] . 

При расчете напряженного состояния представим элемент в виде двух состав- 
ных частей: диска с сердечником (/) и стержневой системы (II) (рис. 7). При этом на- 
пряженное состояние элемента будем рассматривать как сумму плоского осесиммет- 
ричного и линейного напряженных состояний с компонентами а г , ст ѳ , а х . Напряжения 
в диске ( I Определим как функцию давления на поверхности контакта (ж ~ к ). Для 
многослойного диска имеет место зависимость [а ] х [Р ] = [Ь ] , где [а ] - матрица 

податливости, [ Р ] - матрица контактных давлений по границам слоев, [I ] - матрица 
нагрузок [8]. 

Для двухслойного диска имеем 


К]х[Р] = [І], 


( 1 ) 



2 ( „ 2 А „ 2 ( іЛ 



3 


Рис. 6. Структурная модель легированной стали: 

1 - матрица (углеродистая сталь); 2 - карбиды легирующих элементов; 
3 -элемент объема материала 
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Рис. 7. К расчету напряженного состояния элемента материала: 1 - цилиндр (уг- 
леродистая сталь - ж); 2 -сердечник (карбиды - к); /- диск с сердечником; 
//-стержневая система 


Поскольку объем матрицы значительно превосходит объем наполнителя, то при 
расчетах напряжений можно считать, что диск имеет бесконечную толщину, т. е. 

= 2 . Из (1) имеем 


г -> со , а Ііша 
*• 1 1 


г 2 ->оо 


Р\ = 


(“*■ -а )Е Е АТ 


ЯДц„-1)-.ЕЛц,+1)- О) 

Напряжения в диске бесконечной толщины определяются выражением, 


При 


Г = Г] а* ѳ = ± Рі , сг* ѳ = - Р] : 


СТ г,Ѳ — 


а- -± (Ѵж-аЖЕуАТ 
'■* ■е,(р ж+ 1)-Е,(р 1 -1)' 


_ (а ж .-а К )Е К Е Ж АТ 


Осевые напряжения в стержневой системе (II) определим при условии равенства 
осевых перемещений ее элементов М ж = М к , При этом имеем М ж = -ТУ" = // или 
(а к ~а ж )Е к А к АТ 

Для цилиндра бесконечной толщины 


# = 


к ж к к 

1 + АА 


Я -А. 


ІішА'1 ѵ ( а к а .ж ) Е > Л А' Г и < ~~Г~’ г Д е - площади попереч- 
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ных сечений цилиндра и сердечника, соответственно. Окончательно получим 
°:=(а ж -а к )Е к АГ, <=0. 

Сравнение результатов эксперимента и расчета на основе предложенной модели 
показали их удовлетворительную сходимость (Л^ сч =6984 тц, М жсп =5000 тц). При 
этом в качестве материала матрицы принималась сталь 15 (с таким же содержанием уг- 
лерода, как и сталь 12Х18Н10Т). Превышение N ^ над N жсп может быть объяснено 
наличием остаточного ресурса у образцов после термоциклирования. 

Так, испытания образцов при растяжении после 5000 тц показали снижение ис- 
тинного предела прочности стали 12Х18Н10Т на 60-70 %. Проведены сравнительные 
практические расчеты циклической долговечности ствола 30мм артиллерийского авто- 
мата при режимной стрельбе с темпом 600 выстрелов в минуту и длиной очереди п і 
(і— 5,10,15,20 - количество выстрелов в очереди) с использованием обоих вышерас- 
смотренных подходов. Расчеты проводились для наиболее нагруженного сечения ство- 
ла (0,15 м от начала калиберной (нарезной) части канала ствола). Материал ствола - 
сталь 38ХНМФА с категорией прочности КТ 80 (ст О02 =800 МПа). Материал матрицы 
сталь 40 ( а и = 422 МПа). 


Из рис. 1 видно, что интенсивность напряжений в матрице ((Мизес Ж) меньше, 
чем в наполнителе (карбиды хрома, никеля, молибдена — Мизес К) приблизительно в 2 
раза, в то время как характеристики прочности указанных карбидов превышают соот- 
ветствующие характеристики углеродистой стали в 5 — 10 раз. Отмеченный факт ука- 
зывает на необходимость использования при оценке долговечности конструкции харак- 
теристик прочности и пластичности матрицы. 

На рис. 8 представлено распределение интенсивности механических напряжений 
в сечении нарезного ствола ("Мизес показал. " - интенсивность напряжений в % 

° та о,о2 X а на Р ис - 9 - распределение интенсивности термомеханических напряжений 
при макроструктурном подходе. 


Для оценки долговечности ствола воспользуемся уравнением вида [9], [10]: 


е „ = 


4А 


• + 


1 -і 




а 


- г 


т П<1 ЮО 

здесь е, = 0,5 1п , е : 

100- ѵд 

аі, = 0,4ст и ; для сталей с 
сД, = (0,54 - 0,0002)ст„ . 


1 + р 

~^" а > ; для сталей с о и <700 МПа (сталь 40), ш=0,5 и 
ст„> 700 МПа (сталь 38ХНМФА) т=0, 36+0.0002), 
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Рис. 8. Распределение интенсивности механических напряжений в сечении нарез- 
ной части канала ствола 


■А — Х4р«і (2-9) 



Н«м*р* ташкиктр*» 


Рис. 9. Интенсивность напряжений на поверхности нарезной части ствола 


При г = = 0 (от нулевой цикл) имеем е 

8 

шах, 


На рис. 10 представлены результаты расчета долговечности на основе макро- 
структурного и микроструктурного подходов. 
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Рис. 10. Сравнительные результаты расчета долговечности при микроструктурном 
(а, б, в) и макроструктурном (в) подходах 


ВЫВОДЫ 

1 . Показан возможный уровень мартенситного и аустенитного превращения в материа- 
ле орудийной стали 38ХНМФА в условиях стрельбы короткими и длинными очередями 
при темпе стрельбы 600 выстрелов в минуту, а также величины радиального термиче- 
ского износа. 

2. Разработан метод оценки циклической долговечности и механической живучести 
конструкций из высоколегированных сталей, в основу которого положена двух компо- 
нентная структурная модель исходного материала, а также - результаты численного 
анализа, полученные на ее базе, которые показали удовлетворительную сходимость с 
результатами экспериментальных исследований при изотермическом термоциклирова- 
нии стали 12Х18Н10Т. 
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МЕТОДИКА РАСЧЕТА УПРУГИХ ПАРАМЕТРОВ НАНОЭЛЕМЕНТОВ 

А. В. ВАХРУШЕВ, А. А. ШУШКОВ 

Институт прикладной механики УрО РАН, Ижевск, Россия 

АННОТАЦИЯ. Разработана методика расчета упругих констант наноэлементов 

методом согласования перемещений составляющих их атомов с полями перемеще- 
ний упругих “эквивалентных” элементов при различных видах нагружения. Приве- 
дены результаты расчетов модуля упругости наночастиц, нагруженных растяги- 
вающими сосредоточенными силами. Получены расчетные зависимости модуля 
упругости Юнга от диаметра наночастиц. 


ВВЕДЕНИЕ 

Известно, что при изменении характерного размера наноэлементов (нановолокон, на- 
нотрубок, наночастиц и т.д.) их физико-механические характеристики: прочность, мо- 
дуль упругости, деформационные и другие параметры, - изменяются на порядок 
[1,2,4, 5]. Расчетные и экспериментальные исследования показывают, что это обуслов- 
лено, в первую очередь, существенным изменением (причем не обязательно монотон- 
ным) структуры и формы наноэлемента [8-11]. Экспериментальное определение меха- 
нических характеристик наноэлементов представляет технически сложную и трудоем- 
кую задачу, вследствие малости размеров данных объектов. Поэтому моделирование 
является альтернативным и перспективным способом установления данных характери- 
стик наноэлементов [3,6,7]. 

Целью настоящей работы являлась разработка методики расчета упругих констант 

наноэлементоа и анализ их изменения в зависимости от характерного размера нано- 
элемента. 

ПОСТАНОВКА ЗАДАЧИ 

Рассмотрим наноэлемент объёмом О и поверхностью 8, состоящий из N атомов (рис.1 - 
1), и подобный ему по форме и размерам “эквивалентный” упругий элемент (рис.1 -2), 

нагруженные одинаковой системой уравновешенных внешних поверхностных/^ и 
объёмных/^ сил (Рис. 1). 
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Рис. 1 . Наноэлемент (1) и упругое “эквивалентное” тело (2) 


Введём определение “эквивалентного” упругого элемента. Материал данного эле- 
мента тела считаем однородным и изотропным, следовательно, его упругие свойства 
определяются двумя константами [14]. Следует отметить, что в общем случае свойства 
наноэлемента могут быть неоднородными (возможную анизотропию свойств в данной 
работе не рассматриваем). Это обусловлено изменением атомной структуры наноэле- 
мента в области его поверхности [7,9,11]. Поэтому понятие упругого “эквивалентного” 
элемента основывается на одинаковом отклике наноэлемента и упругого элемента на 
заданную систему нагрузок. Иными словами — изменение формы и размеров упругого 
эквивалентного элемента и наноэлемента должны быть равными при действии оди- 
наковой заданной системы нагрузок. 

Деформирование наноэлемента будем рассчитывать методом молекулярной дина- 
мики, согласно которому движение системы атомов, входящих в наноэлемент, опреде- 
ляется системой дифференциальных уравнений 


а 2 г. Л - ~ иг - 

т і , і = ~Ъ Щ ) + /і (0 - ' - р , 

Ш 7=1 Ш 


(1) 


где т і ,Т' І Г- (х ; , у і , 2 { ) - масса и радиус- вектор /-го атома, соответственно; 

г і] ~ ( х і — х .і ) +ІУі~ У]) + ( 2 і~ г у) - радиус - вектор между атомами іи у; 
■^(Д/ ) " сила взаимодействия между атомами іи ^ \ Ъ і - коэффициент “трения”; 
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Р ~ + Р* - суммарный вектор внешних сил; /КО" случайная сила, действующая 

на /-ый атом, определяемая из распределения Гаусса со следующими свойствами: сред- 
нее значение случайной силы 7/(0 равно нулю; она не коррелирует со скоростью 
Ц (0 рассматриваемого атома, так что < У/(0^Л'(0' > равно нулю и 

< 7і{*)/і(®) >= 2 квТ 0 Ьі т і5{і). Здесь к в - постоянная Больцмана; ^(/) - дельта функ- 
ция Дирака; Т 0 - температура. 

Как правило, сила взаимодействия между атомами определяется соотношением 


дѴ(г и ) 

Г(^) = - , (2) 

Ч 

где V (Гу ) - потенциал. 

Постановка задачи статического расчета упругого изотропного тела включает: 

- дифференциальные уравнения равновесия 

V *<т = 0, (3) 

уравнения связи между тензором деформаций 8 и вектором перемещений и 


8 = ^-(и ® V + V ® и), 
определяющие уравнения упругого тела 


(4) 


а = 0(х)е, ( 5 ) 

где Б (х) - матрица упругих констант материала. 

Решение системы уравнений (3)-(5) позволит вычислить вектор перемещений, 
тензор деформаций и тензор напряжений в упругом “эквивалентном” элементе, кото- 
рые зависят от матрицы упругих констант. Согласование решений задачи молекуляр- 
ной динамики и теории упругости осуществляем по векторам перемещений в точках, 
совпадающих с положением атомов нанообъекта 


Ы е 1і тсІ ’ (6) 

где " векторы перемещений упругого “эквивалентного” элемента и наноэле- 

мента, соответственно. 
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Для выполнения условия (6) осуществляется варьирование упругих констант уп- 
руіою эквивалентного элемента, таким образом, чтобы суммарная ошибка, опреде- 
ляемая разностью векторов перемещений, была минимальна. 

МЕТОДЫ РЕШЕНИЯ И РЕЗУЛЬТАТЫ РАСЧЕТОВ 


В качестве примера, поясняющего разработанную методику, приведем расчет модуля 
упругости наночастицы, нагруженной осевыми растягивающими силами (рис. 2). 


В расчетах использован потенциал Леннарда-Джонса 
Щ) = 4б 


( V 2 

а 


\ г ч) 


( Л 6 

(У 


V г а ) 


(7) 


Здесь Б - глубина потенциальной ямы; <т - значение Ту , при котором V ( гг ) равно ну- 
лю. 

Для интегрирования системы уравнений (1) и (4) используем алгоритм Верле 
[12,13]. 

Рассмотрим использование предложенной методики на примере наночастиц цин- 
ка. Количество атомов в наночастицах варьировалось от 216 до 9261. Диаметр равно- 
весных наночастиц цинка при этом принимает значения от 12 до 50 ангстрем. 

Величина силы растяжения составляла =2.086*10 ' ' Н и подбиралась та- 
ким образом, чтобы атом, к которому она прикладывается, не отрывался от частицы. 
Основным результатом расчета деформирования наночастицы является зависимость 
радиальных перемещений атомов и от расстояния до центра масс наночастицы г 
(рис.З). 



Рис. 2. Равновесная наночастица (1) и упругий “эквивалентный” элемент (шар) 
(2), растягиваемые сосредоточенными силами 
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Рис. 3. Зависимость перемещений и от расстояния г до центра масс наночасти- 
цы цинка, содержащей: а - 512, Ь - 3375, с - 8000 атомов, соответственно; 
и иг приведены на графике в ангстремах (А) 


Задача расчета упругого эквивалентного” элемента (шара) при растяжении со- 
средоточенными силами имеет аналитическое решение [14]: 


и = 


(1 + ѵ)Д 


2лг 0 Е 


X (4л + 3)|(соза) Р 2«+і( С08а ) + [„(2л + 1 )р 2 - п(2п - 1) - ѵ(4л + 1)! 81 

п= о I 2 п + 1 


* 3,+д 


(соза) 


т у 


У 


2п+\ 


( 8 ) 


где г 0 - радиус шара; Е - модуль Юнга; ѵ - коэффициент Пуассона; Р 2п ( соза) - по- 
линомы Лежандра; р = г/г 0 < 1 ; Е 1 - величина растягивающей силы. Функция т п оп- 
ределяется из выражения: 


т п =( 1 + ѵ)(4л + \)п + п + \ 


(9) 
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Рис. 4. Зависимость перемещений и от радиуса г: 1 - для упругого шара. 
2 - для наночастицы, состоящей из 2197 атомов;а) Е = 9.5* 1 0 10 Па 
Ьб) Е= 5.25*10'° Па 


На рис.4-а-1 сплошной линией представлена зависимость перемещений от рас- 
стояния до центра "эквивалентного” упругого элемента для модуля Юнга 

Е = 9.5 * ІО 10 Па и коэффициента Пуассона ѵ = 2.1 . На этом же графике (рис.4-а-2) 
приведена зависимость радиальных перемещений атомов и а от расстояния до центра 
масс наночастицы г. Видно, что указанные зависимости не совпадают. Поэтому, изме- 
няя модуль Юнга, добиваемся слияния данных кривых (рис. 46). Критерием этого явля- 
ется среднеквадратичная ошибка, изменение которой в зависимости от модуля Юнга 
приведено на рис. 5. Из данного графика видно, что среднеквадратичная ошибка имеет 
ярко выраженный минимум, которому соответствует наилучшее совпадение векторов 
перемещений. В данной точке и определяется модуль упругости наночастицы. 

Выполняя данную процедуру для наночастиц различного диаметра, строим зави- 
симость модуля Юнга от диаметра наночастицы (рис.5). Очевидно, что модуль Юнга 
увеличивается при уменьшении размера наночастиц, и это качественно совпадает с 
экспериментальными данными [1]. 

По аналогичной методике построены расчетные зависимости модуля Юнга от ра- 
диуса наночастиц цезия и магния. На рис.7 приведена зависимость относительного мо 

дуля Юнга Е от радиуса г для наночастиц из различных материалов. В качестве от- 
носительного модуля выбрано отношение модуля Юнга к асимптотическому его значе- 
нию при максимальном диаметре наночастицы. 

Рис.7 иллюстрирует общую тенденцию увеличения модуля Юнга, однако у раз- 
личных материалов она неодинакова. Относительный модуль цинка увеличивается бо- 
лее чем в 6 раз, а модуль цезия только в 1 .5 раза. 
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0.0 1.0 2.0 3.0 4.0 5.01 6.0 7.0 8.0 9.0 10.0 


Рис.5. Зависимость среднеквадратичной ошиоки \Ѵ от модуля Юнга Е упругого 

эквивалентного элемента; ^ 2 і 97 " модуль Юнга (Па) наночастицы, со- 
держащей 2 1 97 атомов 
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Рис.6. Зависимость модуля Юнга Е (Па) от радиуса г (А) наночастиц цинка 


Е 10' 11 






♦ 




















▼ 






♦ 







1 



♦ — 

^ 






# | 

♦ 





Г 


ХИМИЧЕСКАЯ ФИЗИКА И МЕЗОСКОПИЯ. Том 7, № 3 


283 


А. В. ВАХРУШЕВ, А. А. ШУШКОВ 



Рис. 7. Зависимость относительного модуля Юнга Е от радиуса г (А) для нано- 
частиц из различных материалов: 1- цезий, 2-магний, 3- цинк 


ВЫВОДЫ 

Разработана методика определения модуля упругих констант наноэлементов. Прове- 
денные расчеты позволили выявить основные закономерности зависимости модуля 
Юнга от размера наночастиц. Численные расчеты показали, что модуль Юнга увеличи- 
вается при уменьшении размера наночастиц цинка, цезия, магния. 
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А.А. ВАХРУШЕВ, А.М. ЛИПАНОВ, А. В. ВАХРУШЕВ 
Институт прикладной механики УрО РАН, Ижевск, Россия 


АННОТАЦИЯ. Рассмотрена задача моделирования двумерных течений неньюто- 
новских жидкостей в криволинейных каналах. Проведены численные расчеты па- 
раметров псевдопластической среды, описываемой с помощью реологической мо- 
дели Пауэлла - Эиринга, и исследованы физические явления, возникающие при ло- 
кальном изменении температуры стенки канала. 


ВВЕДЕНИЕ 


В процессе движения жидких сред в различных устройствах, таких как трубопроводы, 
перерабатывающее оборудование, при извлечении вязких ископаемых из несущих пла- 
стов зачастую геометрия каналов, по которым происходит движение, является опреде- 
ляющим фактором, влияющим на основные параметры потока. Наличие выступов, из- 
гибов приводит к возникновению вихревых течений, уменьшению расхода жидкости 
вследствие возрастающего градиента давления. Дополнительными факторами, услож- 
няющими анализ указанных физических процессов, являются реологические характе- 
ристики жидкости: вид реологического закона, определяющий, каким образом реагиру- 
ет среда с той или иной внутренней структурой на внешнюю нагрузку; пределы изме- 
нения коэффициента вязкости, его связь с температурным полем и т.д. Многофактор- 
ное исследование комплексного влияния указанных выше факторов на параметры те- 
чения обусловливает актуальность использование математического моделирования 
астоящее исследование является развитием работы [I]. в которой рассматривалось 
течение неньютоновских жидкостей в прямом канале при наличии структурных (дис- 
кретных) изменений вязких свойств среды под действием температурного поля Рас- 
смотрена задача течения в плоском двумерном канале, состоящем из прямых участков 
соединенных изогнутым ("коленом”), в предположении, что на входе в изогнутую 
часть на одной из адиабатических стенок происходит нарушение теплоизоляции, и 
жидкость локально охлаждается. 
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ПОСТАНОВКА ЗАДАЧИ 

Для расчета параметров течения в комбинированном канале, состоящем из прямоли- 
нейных и криволинейных участков, необходимо использовать криволинейные коорди- 
наты для аппроксимации конечных разностей в криволинейной части трубы. В общем 
случае применяются различные методы построения разностных сеток - аналитические 
или численные [2,3]. Общее требование к различным методам построения - линии сет- 
ки должны быть ортогональными для повышения точности расчетов. В случае сложных 
нестационарных задач требуется перестроение сетки на каждом временном шаге, что 
приводит к увеличению объема временных ресурсов, необходимых для решения зада- 
чи. Предположим, что изогнутая часть канала имеет постоянный радиус кривизны. В 
этом случае можно использовать полярные координаты г,ср. Они являются частным 
случаем цилиндрических координат при отсутствии движения вдоль оси 2 [4]. Исполь- 
зуя зависимости 


X — Г • С05 <р, Г ~ X х + у ’ 
у^г-йпср, ^, =агс1ап >:. 

X 

можно получить аналитические формулы перехода от уравнений Навье -Стокса в де- 
картовых прямоугольных координатах к их записи в криволинейных координатах. 

Общая форма уравнений Навье -Стокса с учетом переменной вязкости может 
быть представлена в виде [7]: 


Дѵ,_ 1 дР | І Э(^) | 1 5(г 2 -^) 
Оі рг д(р г д(р г 2 дг 

Оѵ. _ 1 дР { 1Э(^,) | і а(г.^) 5, 
йі р дг г д(р г дг г 


\_ 

г 



дер 


дг-ѵ г 

дг 


= 0 , 


где - компоненты тензора напряжений. 


(1) 

(2) 

(3) 


$„=217 


г д<р г 
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<рг 


8 ™ = п 


г дер дг 


<р 

г 


8 гг = 27 


дг 


( 4 ) 


Если предположить, что вязкость жидкости т] не меняется но времени и внутри 
области интегрирования, то уравнения (1),(2) упрощаются до вида: 


дѵ т 5ѵ„ ѵ дѵ ѵ ѵ 

?_4- Ѵ | г <Р 


ді 


дг г д(р 


1 дР 
рг дер 


+ Т] 


Л ѵ <р 2 <Эѵ 
Дѵ — г-н - 

ѵ г 1 г 1 д(р 


( 5 ) 


где 


дѵ, ѵ а дѵ.. ѵ, 


дѵ 

— - + ѵ г —^ + 
ді дг г дер 


<Р г У 


\_ дР 
р дг 


+ П 


Аѵ, 


г 2 дер 


( 6 ) 


л 

А = — г + 


+ - 


2 ' 


дг 1 г дг г 2 дер 

Переход к безразмерным переменным при решении задачи течения 
криволинейном канале позволяет записать (5)-(6) 


( 7 ) 


жидкости в 


_ дР 


+ ■ 


Ж дер Ке 


Д Ѵ + 

<Р 2 2 п 

г г дер 


( 8 ) 


±\__дР 

Ж дг Ре 


ѵ 2 дѵ 

АѴ Г- 

г г д(р 


( 9 ) 


где Ке 


. рУ а Щ-В і) 


- число Рейнольдса, характеризующее тип течения (ламинар- 
ное/турбулентное); К } ,К 2 - внутренний и внешний радиусы кривизны канала; V, - ха- 
рактерная скорость потока, втекающего в канал (в частности, берется максимальное 
значение скорости на оси канала, соответствующее профилю Пуазейля); т), р - кине- 
матическая вязкость жидкости и плотность среды, соответственно. 

Уравнение неразрывности (3) для использования в дальнейших преобразованиях 
можно представить как: 


дѵ <Эѵ 

— - + ѵ+г — - = О 
дер дг 


(Ю) 
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с учетом условия г * О . 

Уравнения энергии движения (8)-(10) записаны для физических переменных 
скорость-давление’ . В первом приближении будем считать, что расход жидкости на 
входе в канал постоянный. Тогда из расчетов можно исключить давление и свести ис- 
ходную систему уравнений к двум неизвестным - ”вихрь-функция тока”. 

Для вывода уравнения переноса “вихря” со использовались соотношения: 


со 


поляр. 


1 

у 

'«У 

1 

дѵ г 

г 

дг 

дер 


( 11 ) 


ду'поляр ѵ 1 д^поляр. 

дг г дер 


( 12 ) 


со. 


поляр 


ЗѴполяр. , 1 ду/ 1 Э> 


(13) 


Легко показать, что переменные ѵ (р ,ѵ г , о поляр , ,у/ поляр в пределе, при радиусе кри- 
визны г —> оо , совпадают со своими аналогами в декартовых координатах и, ѵ,со,ц/ . 
Это утверждение помогает осуществить переход на границах прямого и изогнутого 
участков трубы от полярных к декартовым координатам. Далее опустим нижний ин- 
декс и будем обозначать “вихрь” и “функцию тока” как со и , соответственно. 

Для получения уравнения переноса вихря в полярных координатах исключим 
давление из уравнений закона сохранения импульса (5)-(6)., используя представление 
со в виде (11): 
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г дер дг 
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г дердг 
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(14) 


Для проведений численных расчетов важной является недивергентная форма 
уравнения переноса. Она позволяет использовать метод характеристик и численные 
схемы, использованные автором при расчете параметров течения в прямом канале [1]. 
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Было получено уравнение переноса “вихря” со в недивергентной форме: 

Рсо _ дсо 1^( ѵ р '<у) 3(ѵ г -<д) 

Рі 3/ г дер дг ( 1 


Полученное уравнение (15) не учитывает диссипативные силы, связанные, в 
случае течения вязко -пластической жидкости, с силами вязкого трения. Для получения 
Источниковых членов в нем были использованы уравнения (1)-(2), в предположении, 
что вязкость жидкости Г) не меняется по пространству и времени. Это приводит (для 
безразмерных координат) к формуле 


Ф 


( о 


Ке 


д 2 со 1 дсо 

7" Н К 

дг~ г дг 


1 д 2 со 
г 2 дер 2 


- Аео. 


(16) 


Выражения для тензора напряжений (4) в цилиндрических координатах исполь- 
зовались для получения диссипативной части уравнения переноса “вихря” в виде: 

Ф« =Д(*7-еО + ф.. 

Такая форма диссипативного члена обусловлена тем, что с помощью нее можно 
построить численную схему, в которой большая часть составляющих вычисляется на 
новом временном слое. В ходе численных расчетов было выяснено, что при этом реше- 
ние обладает большей устойчивостью, что позволяет брать больший шаг по времени; 
скорость расчета возрастает. 

Таким образом, получена диссипативная составляющая уравнения переноса в 
цилиндрических координатах: 


Ф® =А(т]-со)-~- д V 


2 ЗЭ; 
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д 2 г] л 


1 ді) 
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(17) 


Необходимо отметить, что знаки перед слагаемыми в формуле (17) могут ме- 
няться, в зависимости от предполагаемого положительного направления 
“завихренности” со. 

Для расчета параметров течения в изогнутом канале, с учетом реологической 
модели Пауэлла - Эйринга, необходимо записать уравнение энергии в полярных коор- 
динатах. Воспользуемся общей формулировкой уравнения энергии для полярных коор- 
динат [7]: у 
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11 дТ 1 дТ 

где ^ = , = компоненты вектора плотности а энергии (в со- 

Ре г д(р Ре дг 

ответствии с законом Фурье). 

С учетом соотношений для тензора напряжений (4) и переходя к недивергентной 
форме записи уравнения энергии, получим 


дТ | 1 д(ѵ ѵ -Т) 
ді г д(р 
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д(ѵ г -Т) ^ѵ г -Т 1 
дг г Ре 


АТ + 7 • Ф 7 . 


(19) 


где диссипативная функция, определяющая рассеяние механической энергии движения 
жидкости и переход ее в тепловую, вычисляется как 


ф Г = 2 


дер 


■ + ѵ. 
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г — 
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'г 
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( 20 ) 


Очевидно, что для переменной и постоянной вязкости уравнение энергии имеет 
один и тот же вид (в предположении, что остальные свойства среды постоянны по всей 
области интегрирования). 


РЕЗУЛЬТАТЫ РАСЧЕТОВ 

В качестве расчетной области рассматривалась изогнутую часть канала (рис.1), которая 
на практике может представлять собой колено трубопровода, соединяющего два участ- 
ка трубы. Наиболее интересный случай представляет вариант, когда происходит охла- 
ждение жидкости на искривленном участке канала. При этом, помимо вихревых обра- 
зований в течении, снижающих скоростные характеристики потока, возникают допол- 
нительные потери энергии при транспортировке жидкости из-за охлаждения жидкости 
и повышения вязкостных характеристик среды. 
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Рис 1 . Расчетная сетка для моделирования течения в изогнутом канале 

При тестовых численных расчетов было выявлено, что границы изогнутой части 
канала являются источником интенсивных возмущений, поэтому расчет параметров 

движущейся среды целесообразно проводить с присоединением к изгибу прямолиней- 
ных участков (рис. 2), 

Прямые участки канала являются продолжением искривленной области интег- 
рирования с обоих торцов. Как видно из иллюстрации, в этом случае границы расчет- 
ной области достаточно удалены от колена, и возмущения, возникающие в потоке при 
его прохождении, слабо сказываются на входном и выходном участках трубы. 

Конечно-разностный аналог уравнения переноса “вихря” а запишется следую- 
щим образом: 



_ 1 _ 

Ке 




( 21 ) 



Рис.2. Расчетная сетка для моделирования течения в изогнутом канале 
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где 
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Ѵй) = 

* ч г } Аср 


[ а и ■ т >, + а и ■ ~ а о ' Ѣ-м ~ < ' % 


(22) 


г Н = + ^ ! ®«.' -д,” ■«>,]. 


(23) 

(24) 


2,з (л) і.і+\ 03 и -\ 


Цсо\ = 

чі 


г Аг 


+ ■ 


Дг 


2 I *"//+] ^ ' Щ] "*■ Щ )- 1 ] • 




(25) 


В представленных зависимостях (22) и (23) величины <Хц и - характеристики 
компонент скорости (ѵ р и ѵ г , соответственно) через грани контрольного объема [5,6]. С 

у четом данных характеристик, частные производные по пространству, стоящие в кон- 
вективной части уравнения переноса, вычисляются с учетом направления движения 
жидкости. Это делает смысл конечно-разностного уравнения (21) более физическим, по 
сравнению с использованием центральных разностей при вычислении частных произ- 


водных 


Ѳѵ.-Р 


И 


д<р “ дг 

Для решения системы линейных алгебраических уравнений, получаемой при 
использовании численной схемы (21), воспользуемся процедурой АЭІ (Акетайѵе Эі- 
гесбопз Ітріісіі) [2]. В соответствии с данным методом, разобьем задачу на три времен- 
ных слоя: п, п + \ , п + 1 , - а также разделим тангенциальное и радиальное направления: 


. , т 1 І" + ! Д/ , ■> |Н+І 

-+д/гД 


(26) 


П ~ 1 , л . , | | /;+ І Аі г 2 |Я+І п+ і 


со'Г + АН' соГ- — 1 2 


н+І 

СО\ - со 

" Не г 1 и " 


(27) 


В результате получим две трехдиагональные системы уравнений вида 
а ?-а>і+ѵ+ЬГ-а> 0 +с?-щ_и=/?, 


(28) 


а г го>,иі +Ь г ,-со„ + 


-./ г и і ' Ш ІІ -и с ) ' Ш и-\ = //, 


(29) 


где 
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К = 
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1 А I 
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(30) 


коэффициенты трехдиагональной матрицы при разложении уравнения переноса со 
вдоль тангенциальной направляющей ср ; 


а г 

- — д 12 

1 А 1 

1 А/ 

,і 

Аг ■' 
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= я 2 ' 

1 А( 

1 А 1 


А г ' Ке Аі'~ 2-Кег^’ 

— и+А 

іі ~ (о ц \ 


(31) 


коэффициенты трехдиагональной матрицы при разложении уравнения переноса со 
вдоль радиальной координаты направляющей г . 

Полученные трехдиагональные системы конечно-разностных уравнений 
(28), (29) решаются методом “прогонки” [5,6], описанным ранее. 

Уравнение Пуассона (13), связывающее “завихренность” со с “функцией тока” 
в полярных координатах, можно представить в конечно-разностном виде, используя 

центральные разности. Данные аппроксимации частных производных дают второй по- 
рядок точности по пространству. 


г і Аг 




[^ + и-2-^+^ М/ ] = 


со іг 


(32) 


Решение системы линейных алгебраических уравнений выполнено методом по- 
следовательной верхней релаксации [2,5,6]. Для этого, предварительно сгруппировав 
соответствующие слагаемые, приведем равенство (32) к виду 


Ги = 


2 - о »; +0 


■ [умі + Ѵі-Уі + (Д; + «,-)■ у/ и+] + (Д 2 




( 33 ) 
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где коэффициенты а,, В, вычисляются как а = — г. -Дю, В - — - ^ . 

7 7 2 7 7 Дг 

Сложим правую часть равенства (33) с тождеством — ц/ ^ = 0 и перегруппиру- 

ем слагаемые 


Г» + 2 (/? 2 + 1 ) ' + ^ + +а ' У Ѵи*' + 

+ (Р] - «/ ) ■ ѴѴі - 4 «>,у - 2(^- + 1) • УУ, 


(34) 


В итоге, предполагая, что в ячейках с индексами /-Іи у 
ны на новом п + 1 временном слое, получим 


- 1 значения извест- 


V 


= Гч + 


' \у'і + \,і + Ѵ'му + (Р) + а і ) ■ ^,"+1 + 


•(А 2 + 1 ) 

+ (Д! -а, - 4 а)ар - 2 (^+ 1 )^”], 


(35) 


где г - релаксационный параметр, для которого выполняется условие: 0 < г < 2 . 

Неявная схема обладает лучшей сходимостью при решении уравнений, описы- 
вающих движение сред со сложной реологией (рис.З). 

Если использовать численную схему для недивергентной формы уравнения пе- 
реноса вихря, то, как показали численные расчеты, погрешность такой схемы достаточ- 
но велика, а ошибка вычислений быстро накапливается. Схема обладает слабой устой- 
чивостью, что требует уменьшения используемого шага по времени, а 



Рис 3. Реологическая поверхность, описывающая поведение псевдопластической 
жидкости 
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полученные результаты не соответствуют физической природе процессов. Так, напри- 
мер, в области повышенной вязкости возникает источник движения, решение теряет 
устойчивость, что приводит к “аварийному” завершению расчетов. 

Для расчета течения по составной трубе потребовалась стыковка уравнения пе- 
реноса вихря на границе участков. Сравним формулы для “завихренности” со в де- 
картовых и полярных координатах: 

вѵ ах' дг г г а ѵ 

По размерности обе величины совпадают, а с учетом того, что указанные части 
канала стыкуются под прямым углом, на граничной области = и, ѵ г = ѵ . Также, при 
выборе достаточно маленьких пространственных шагов расчетной сетки Дх, Д у, АсриАг 
в пределе будут выполняться соотношения г-А(р-+Ах, Ау = Аг. Тогда очевидно, что 
основной вклад в выражении Асо = со поляр - « декарт , основываясь на (36), будет состав- 
лять Ѵр/г . Введем обозначение новой функции со 


= ^поляр ~б, где д = — , 


(37) 


которой соответствует новое уравнение переноса 


дсо 1 дѵ ѵ -а дѵ г -со 1 д 2 т]-со 1 дт]-сд д 2 п-сд 

ді г д(р дг г 2 д + г д + + дг 2 (^) 


где диссипативная функция вычисляется с помощью соотношения (17), а С- опре- 
деляется как 


дб 1 дѵ -8 дѵ г -5 
ді г д(р дг 


1 д 2 г/-3 1 дг)-8 

— — + — + 

{г 2 д г д 




(39) 


С помощью формул (37)-(39) проведена серия расчетов, как с учетом постоян- 
ной вязкости среды (рис.5), так и при наличии зоны охлаждения и повышении вязких 
свойств жидкости (рис.6). 

На основе экспериментальных данных известно [8], что при течении вязкой 
жидкости по составному каналу (рис. 2), возникают вихревые течения. Их форма, на- 
правление движения и области образования отображены на рис. 4. 
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Рис. 4. Возникновение вихревых течений при движении жидкости по изогнутому 
каналу 



с 


й 


Рис. 5. Течение в изогнутом канале при постоянной вязкости 
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Рис. 6. 1 ечение в изогнутом канале при наличии скачка вязкости 

Данное явление наблюдается и в численном эксперименте (рис.5): обратное те- 
чение возникает на входе в изогнутую часть трубы и на некотором расстоянии от выхо- 
да жидкости из данного участка. Области обведены на рисунке. 

Для расчета течения вязкой жидкости в изогнутом канале границы расчетной 
области принимались достаточно удаленными (порядка 5-6 поперечных размеров кана- 
ла) от изогнутой части. Это условие использовалось для того, чтобы вихревые зоны. 

образующиеся в потоке, не оказывали существенного влияния на выходные границы 
канала. 

Расчеты проводились так же, как и в задаче течения жидкости по прямому кана- 
лу, при числе Рейнольдса Ке = 10 в невозмущенном потоке, а в локальных областях 
вязкость г] увеличивалась до 50 раз. 

В случае постоянной вязкости (рис.5) поток жидкости, входящий по прямому 
участку канала в изогнутый, “прилипает” к нижней стенке. В результате, в верхней час- 
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ти образуется разряженная зона с вихревым течением. Далее, под действием центро- 
стремительной силы, поток жидкости отклоняется к верхней стенке канала, и на выходе 
из изогнутой части образуется “теневая зона", вблизи внутренней твердой поверхности 
с циркуляционным течением. 

Наибольшие энергетические потери, связанные с вязким трением, наблюдаются 
в потоке в областях, отмеченных сплошными линиями на графике распределения 
“завихренности” со (рис.5(с)). 

Течения с локальными областями изменения вязкости в изогнутом канале обла- 
дают еще более сложной структурой (см. рис. 6). На графике (Ь) данного рисунка пока- 
зано распределение вязкости в случае охлаждения участка нижней твердой границы на 
входе в “колено”. Встречая вязкую “преграду”, поток отклоняется к верней стенке. 
Также как и в случае с течением в прямом канале, наблюдается явление “биения 
струи’ . Эффект данного явления усиливается появлением третьей вихревой зоны (на 
выходе из изогнутого участка). Она представляет собой развившуюся циркуляционную 
область при течении жидкости с постоянной вязкостью за счет дополнительных 
колебательных’ движений основного ядра потока. 

Видно, что силы вязкого трения достигают максимального значения на границах 
области повышенной вязкости, что отображено на рис. 6(с). 

При возникновении скачка вязкости в изогнутом канале наблюдается большее 
число областей с вихревым движением (по сравнению с потоком, проходящим через 
прямой канал). В случае течения с постоянной вязкостью образуется две циркуляцион- 
ные области, а при локальном увеличении вязкости их число увеличивается до трех. 
Эю явление связано с тем, что помимо сил вязкого трения, создающих дополнительное 
сопротивление в области охлаждения жидкости, поток, устремляющийся к стенке кана- 
ла под действием центростремительной силы, “отражается” от твердой границы. В ре- 
зультате образуется “зона подсоса”, обусловливающая появление вихря. 

Образование дополнительных вторичных и третичных циркуляционных течений 
при движении вязкой жидкости по трубам сложной формы приводит к дополнитель- 
ным энергетическим потерям при транспортировке высоковязкой среды: возрастает 
градиент давления вдоль канала; как следствие - снижается расход жидкости. На осно- 
ве построенной методики можно исследовать, как повышение вязкости среды повлияет 
на картину течения и как, например, за счет установки дополнительных нагреватель- 
ных элементов, можно бороться с вихревыми образованиями в потоке. 
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О КРИТЕРИИ БЕЗОПАСНОГО СЖИГАНИЯ ДЕФЕКТНЫХ ЗАРЯДОВ 
ТВЕРДОГО ТОПЛИВА ПРИ ПОНИЖЕННЫХ ДАВЛЕНИЯХ 

Д. А. ПОСПЕЛОВ, В.Ю.МЕЛЕШКО, П.Н.НАУМОВ, Г.Я.ПАВЛОВЕЦ 

Военная академия РВСН имени Петра Великого, Москва, Россия 


АННОТАЦИЯ. Рассмотрены дефекты твердотопливных зарядов снимаемых с воо- 
ружения ракет, влияющие на процесс сжигания. Предложены критерии, позволяю- 
щие оценить влияние дефектов заряда на увеличение поверхности горения и расхода 
продуктов сгорания. Определены условия, при которых растрескивание внутренней 
поверхности заряда можно считать не опасным для его утилизации. 


К настоящему времени можно считать очевидным, что утилизацию зарядов твердого 
ракетного топлива приходится осуществлять при условии наличия в зарядах дефектов в 
виде трещин и отслоений, локализованных в различных местах конструкции. При этом 
накопленный опыт сжигания таких зарядов без соплового блока (давление в камере, 
как правило, не превышает ІМПа) на специальном стенде показал следующее: 

чисто механического опасного развития дефектов в процессе горения при давлении 
в камере сгорания порядка 1 МПа не наблюдалось; 

поверхности горения с обнаруженными отслоениями заряда от корпуса могут быть 
флегматизированы - пропитаны специальными составами, предотвращающими вос- 
пламенение дополнительной свободной поверхности при горении заряда в процессе 
утилизации, механического разрушения зарядов от действия внутреннего давления 
при этом также не происходило; 

при горении заряда могут возникать скачки расхода продуктов сгорания, что явля- 
ется опасным как с точки зрения сохранения работоспособности самого стенда, так 
и соответствия работы стенда экологическим требованиям. Такого типа явления мо- 
гут быть вызваны образованием дополнительной поверхности горения вследствие 
растрескивания топлива при действии на заряд внутреннего давления и воспламе- 
нения топлива по внутренней поверхности образованных полостей. 

Последнее наблюдение представляется достаточно опасным и делает необходи- 
мым выбор приемлемого исходного критерия безопасной утилизации дефектного заря- 
да методом бессоплового сжигания. Такой критерий может быть выбран на основе ог- 
раничения величины образуемой при горении топлива дополнительной поверхности 
горения: 
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где 8 - площадь поверхности горения; 

Д8 - дополнительная поверхность горения на канале заряда за счет проникновения 
пламени в полости дефектов [1]. 

При действии внутреннего давления на жестко скрепленный с корпусом заряд ос- 
новную нагрузку принимает на себя корпус (подкрепляющее действие заряда составля- 
ет величины порядка нескольких процентов), поэтому в первом приближении можно 
считать поле деформаций в заряде заданным. Считая деформации заряда малыми, поле 
деформаций плоским, а материал несжимаемым из решения задачи Ляме получим за- 
висимость тангенциальной деформации в от радиуса г: 


(Г) 


г 


где е к - тангенциальная деформация корпуса; 
а - радиус корпуса. 


При длительном хранении заряда твердое топливо существенно и неоднородно 
охрупчивается, причем деформативность топлива е пра) тем меньше, чем ближе выбран- 
ная область к поверхности канала. Примем, что при действии внутреннего давления 
(т.е. при задании деформации корпуса) радиальные поверхностные трещины достигнут 


вершинами некоторого радиуса г 0 , такого, что 


е М = е п/юі (г 0 ). 


( 2 ) 


Условия (1) и (2) задают некоторое алгебраическое уравнение для определения ве 


личины к : 


а 


8 , 


Если задать деформативность в вершине трещины получим: 



В этом случае при Ъ > г 0 величина глубины трещин / 


примет значение 
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ѵ — Ь ~ Ъ 

г 


\ 

а 1 

'о и и 

А 

^ прс< ) 

Ь 

У 


где Ь - внутренний радиус канала. 

Последнее соотношение представляет собой уравнение относительно неизвестной 
величины - глубины трещин І тр , поскольку деформативность топлива е пра) есть пере- 


менная величина. е прей - е прей (г 0 ), т.е. е пред = е п/К() (і тр ) . Рассмотрим два предельных 

случая: когда трещина расположена вблизи поверхности и величина е прсд почти посто- 
янна, и когда глубина трещины совпадает с глубиной приповерхностного охрупчивания 
топлива. 


Первый случай (трещина расположена вблизи поверхности) 

Величина раскрытия трещин в основании (на внутренней поверхности канала) может 
быть оценена, если известно количество трещин. Пусть N — количество трещин, 
8 пред {Ь) - деформативность топлива на внутреннем канале (примем, что в области рас- 
трескивания заряда она известна), тогда величина раскрытия 5 определяется из соот- 
ношения 


8 = 


2 пЪ 


N 


а 


е ^-е пр А Ь ) 


( 3 ) 


Таким образом, из соотношений (1;2;3) можно определить отношение параметра 
— , который является основным при оценке как возможности воспламенения внутрен- 
ней поверхности дефекта, так и лавинообразного роста давления в полости дефекта при 
условии горения топлива по внутренней поверхности дефекта [1 ;2]. 

Количество трещин можно оценить из энергетических соображений. Действитель- 
но, пусть материал в приповерхностном слое толщиной, равной протяженности трещи- 
ны, разгрузился с выделением некоторой энергии АЭ, которая была затрачена на обра- 
зование новой свободной поверхности ДА т/ , [3]. Пусть у - удельная поверхностная 
энергия, тогда 

(4) 


Зная общую площадь и глубину трещин легко оценить величину N - общее коли- 
чество образовавшихся трещин. 

Для оценки высвободившейся энергии АЭ примем, что в приповерхностном слое 
по внутреннему каналу заряда на глубину трещины деформация топлива после образо- 
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вания трещин остается на уровне предельного значения е п[т) . Тогда, в первом прибли 

женин, без учета трехмерности напряженно-деформированного состояния топлива ве- 
личину энергии АЭ можно оценить из следующих соотношений: 


ДЭ = ЩЕ^Ѵ-Е^2*Ы т ,,1, 


где I - длина заряда, Е - модуль Юнга. 
С учетом зависимости ( 1 ), 


Ш Е У с * Ѵ = 2/гіе [ ~гг(іг = 2лЬЕ^ 
1 2г 


' зоря (к) 


к ЕЕ' 


еУ_\_ 

2 г 2 2 


/>-/, 


'тр 


е 2 а А 

( 

1 

> 

1 

2 

Т 

'--О 

ь 2 ' 

у 


лЬЕв 2 а^ а 4 

— ъ ; І тр Ь = ^ЕеХ р Ь~. 


Таким образом 


АЭ - пЫ^ЬЕ 


С 4 

а 2 2 

Р — Р ^ 

^4 °а- °прс() 


(5) 


Пусть к Е - значение исходного коэффициента запаса заряда по деформативности. 
р 0 - расчетное значение внутреннего давления на канале при работе двигателя, р у - 

максимальное давление в камере сгорания в процессе утилизации, а У пра) - деформа- 
тивность «свежего» топлива, тогда из формул (4;5) следует 


Д-У. _ ^ ЕЫ тг 


/ 2 лЬЕ 


( е ° ѵ Ѵ 

° пред Ру 


К Ро 




пред 


ЕІ р° 

— т Р^ пред 


2 У 


''пред Ру ^ 


V/ 


V к е е пра> Ро 


Ру 

ѴЩТоУ 


/ 

\ 2 Г А 

г 

° пред 


2\ 


Р 4 
\ пред у 


( 6 ) 


Из полученной формулы видна асимптотика роста относительной величины до- 
полнительной площади горения - величины ^ с уменьшением относительной де- 

формативности топлива - величины щрА 

^ пред 
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Формула (6) позволяет оценить предельные значения глубин трещин / , при ко- 
торых уровни увеличения расхода допустимы, т.е. стенд «справится» с дополнитель- 
ными выбросами продуктов сгорания. 

Пусть х ~ предельное значение относительного приращения расхода продуктов 
сгорания. Тогда 


следовательно, 


ДА 


тр 


8 




1 < 

тр — 


2 УХ 


Ее 


пред 


V,) Ру_ 
к(р„ рс<) Р 0 


V 


Г „ X 


ѵЛд»у 


2 / 


ч 2 Л ■ 


' пред 


\ пред 


С другой стороны, величину свободной поверхности А8 тр можно оценить из соот- 
ношения 


А8 тр =2КЫ тр . 


Полученное соотношение с учетом формул (3) и (6) позволяет оценить степень 
раскрытия трещин в вершине 


2/ 


8 *> 2- 


( а 1_ 

^2 ^пры) 


4 у 


( 4 

Я 2 

\ 

2 


Га 2 

\ 

— тЕ 

и 4 * 

- 8 . 

Л/К*) 

У 

Я 

1 04 
і <5 

^ Е пред 

У 


(7) 


Второй случай (глубина трещины совпадает с глубиной приповерхностного 
охрупчивания топлива) 


Пусть И - глубина охрупчивания топлива, І тр = И . Тогда 


~ 2 пЬ а 2 

о = е, — 

N к Ъ 1 


(ь+и) 


4 пЪ а 1 
»г Е К ТТ • 

N Ь~ 


( 8 ) 


В первом приближении 


А Э = 


Ее. 2 а 4 


2 V 


(6 + Д) 4 


2яЫН » 2Ее 2 ~-—2яЫИ. 

к ь 4 ь 
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Аналогично первому случаю оценивается абсолютное приращение площади 


А? = 


2 Ее 2 а 4 И 
У Ь 4 Ь 


2 яЫИ 


и относительное приращение площади 


А? _ 2 ЕЬе К 2 а 4 И 2 
5 у Ъ 4 Ъ 1 


Последнее соотношение позволяет оценить асимптотику роста дополнительной 
площади горения с ростом глубины охрупчивания топлива. К сожалению, получение 
априорных прогнозов по значениям физико-механических характеристик «свежих» ма- 
териалов или искусственно состаренных достаточно проблематично. Есть основания 
полагать, что и величина удельной поверхностной энергии у также существенно 
уменьшается в процессе охрупчивания топлива, как и снижается деформативность. Ес- 
ли учесть, что модуль Юнга может расти, то легко заметить, что все указанные процес- 
сы приводят к росту оценки дополнительной площади горения А8 . 

Зная величину приращения площади за счет трещин, можно оценить количество 
трещин 


2 Ее 2 а 4 И 


N 


у Ъ Ъ 


— 2 пЪЬИ 


НИ 


2пЕИе 2 а 4 
~у V* 


Зная количество трещин, из соотношения (8) можно оценить величину раскрытия 
трещины в вершине: 


4 яИ^а 1 _ 4яНе *%Г 
N е,; Ь 2 ~ 2 кЕИе 2 а 4 
у Ъ 4 


2 у Ъ 2 
Ее к а 2 


(9) 


Оценка (9) является нижней, а оценка (7) - верхней. 

Таким образом, полученные соотношения позволяют оценить уровни площади до- 
полнительной свободной поверхности, образованной за счет растрескивания топлива 
по поверхности внутреннего канала заряда. При выполнении известных условий [1;2]. 
эта поверхность станет дополнительной поверхностью горения, что увеличит расход 

продуктов сгорания. Зная длину трещины и отношение ^ , можно определить сле- 

8 

дующее: 
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1) Произойдет ли воспламенение внутренней поверхности трещин до выгорания 
топлива на глубину трещины, т.е. следует ли считать увеличение поверхности горения 
на величину А5 тр , для чего оценивается <5 тш - величина минимального раскрытия де- 
фектов, при которых достигается на поверхности в полости дефектов температура вос- 
пламенения 


3 ^ = 8 /. 


тр 



Т*-Т н 1 

1_ [К Ро< 
К 


где р т - плотность топлива; 
с - теплоемкость топлива; 

Я - теплопроводность топлива; 
а - коэффициент температуропроводности; 
Т н - начальная температура заряда; 

Т в -температура воспламенения; 


К 

— - газовая постоянная смеси продуктов сгорания; 


и - скорость горения топлива; 

ѵ 0 - средняя скорость молекул в камере сгорания; 

р 0 - давление в камере сгорания [3]. 

2) Произойдет ли лавинообразное развитие дефекта вследствие роста давления в 

полости трещины [1], для чего оценивается максимальный уровень относительного 

г \ 


давления 


Р 


КРо) 


шах 


в полости дефекта: 


'Р_] 

УРо) тіХ р с х р о’ 

где 5 - площадь внутренней поверхности трещины; 

Р - площадь поперечного сечения трещины в месте ее выхода на воспламененную 
поверхность, 

с, - числовой коэффициент [3]. 

В том случае, если 

а) воспламенение внутренней поверхности трещин до выгорания топлива на глу- 
бину трещины не происходит [2], 
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Ь) максимальный уровень относительного давления 


( \ 

Р_ 

Роу 


в полости дефекта не 


превышает предельно допустимого значения [1], 

с) величина образуемой при горении топлива дополнительной поверхности горе- 
М 

ния: ^ не превышает допускаемый характеристиками стенда уровень относительного 
превышения расхода, 

то растрескивание внутренней поверхности заряда можно считать не опасным для его 
утилизации. 
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МОДЕЛИРОВАНИЕ ПРОЦЕССОВ ПОГЛОЩЕНИЯ И ВЫДЕЛЕНИЯ 
ВОДОРОДА ФУЛЛЕРЕНАМИ И ФУЛЛЕРИТАМИ 

А.М. ЛИПАНОВ, А. В. ВАХРУШЕВ, М.В. СУЕТИН, А.А. ВАХРУШЕВ 

Институт прикладной механики УрО РАН, Ижевск, Россия 

АННОТАЦИЯ. Методом молекулярной динамики выполнены расчеты адсорбции, 

хранения и десорбции водорода фуллеренами С 20 , С 60 , С 80 , Сш, С 24 о и С 54 о и сфор- 
мированными из них фуллеритами при различных внешних условиях. Рассчитаны 
кинетика адсорбции молекул водорода фуллеренами; определены термодинамиче- 
ские параметры, необходимые для устойчивого хранения водорода фуллеренами; 
построена структура, образованная молекулами водорода на фуллеренах при хране- 
нии; определены термодинамические условия выделения водорода. 


ВВЕДЕНИЕ 

Поиски альтернативных источников энергии в связи с истощением природных запасов 
углеродных топлив и ухудшающейся экологией весьма актуальны в настоящее время. 
Водород, в этом смысле, является новым перспективным экологически чистым топли- 
вом. Для накопления водорода и его хранения используют различные устройства (ак- 
кумуляторы водорода) конструкции которых, несмотря на постоянное совершенствова- 
ние, не позволяют достичь высокого весового содержания водорода. У промышленных 
аккумуляторов весовое содержание водорода не превышает 5-6 % [1]. В связи с этим, 
следует отметить, что для накопления водорода можно использовать аккумуляторы, 
созданные на основе наноструктурных элементов: фуллеренов, наночастиц, нанотрубок 
и т.д., - обладающих уникальными свойствами по удержанию различных газов [2,3]. В 
настоящее время в этом направлении проводится большое количество эксперименталь- 
ных работ, но нет однозначности результатов вследствие применения различных мето- 
дик проведения опытов. Кроме того, в большинстве случаев невозможно получить 
идентичные образцы наноструктур для исследования. Поэтому использование методов 
математического моделирования для проектирования и оптимизации структуры нано- 
аккумуляторов водорода весьма актуально и перспективно [4-7]. 

В данной работе в развитие работ авторов [8-10] представлено моделирование 
методами молекулярной динамики процессов поглощения, хранения и выделения водо- 
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рода при различных внешних условиях фуллеренами С 20 , С 60 , С 80 , С 180 , С 240 и С 540 и 
сформированными из них фуллеритами. Рассчитаны кинетика процессов адсорбции 
молекул водорода фуллереном; определены термодинамические параметры, необхо- 
димые для устойчивого хранения водорода фуллереном; построена структура, образо- 
ванная молекулами водорода вне фуллерена при хранении и рассчитаны термодинами- 
ческие условия и кинетика выделения водорода, 

ПОСТАНОВКА ЗАДАЧИ 

Для моделирования взаимодействия системы фуллеренов (фуллеритов) с водородом 
рассмотрим расчётную ячейку с периодическими граничными условиями (рис. 1). 

Движение атомов данной наноструктуры определяется системой дифференци- 
альных уравнений 



( 1 ) 



2 


3 


Рис. 1 . Расчетная схема: 1 - фуллерены, 2 - молекулы водорода, 3 - ячейка с пе- 
риодическими граничными условиями 
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при начальных условиях 



( 2 ) 


где п - число атомов, составляющих наносистему; т г масса і-го атома; Г) 0 ,Г; - на- 
чальные и текущие радиус - вектора і-го атома, соответственно; Р| (г) - суммарная си- 
ла, действующая на і-ый атом со стороны других атомов; О - область, занимаемая на- 
носистемой (3); Ѵ і0 , Ѵ і - начальная и текущая скорости і-го атома, соответственно. 

На границах расчетной области (3), (рис. 1) задаются периодические граничные 
условия, согласно [11]. 

Р і (г) определяется соотношением: 



( 3 ) 


где г - {г 1 ,г 2 >--,г п }; Е(г) - потенциал, зависящий от взаимного расположения всех 
атомов, описывающий связанные и несвязанные взаимодействия. 



( 4 ) 


где слагаемые отвечают следующим типам взаимодействий: Е^ - химическим связям, 
Е ѵ - валентным углам, Е ф - торсионным углам, Е^ - плоским группам иЕ и - Ван - 
дер - Ваальсовым контактам. 

Валентные длины поддерживаются за счёт потенциала: 


Е ь = ^КДЬ-Ь 0 ) 4 , 


( 5 ) 


Ьопсіз 


где К г - константа растяжения связи; Ь 0 - равновесная длина связи; Ь - текущая дли- 
на связи. 

Валентные углы задаются потенциалами 


Е ѵ = ^К ѳ (Ѳ-Ѳ 0 ) 4 , 


( 6 ) 


ап§1ез 


где К ѳ - силовая константа изгиба; Ѳо - равновесная величина угла связи; 
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Ѳ - текущая величина угла связи. 

Энергия торсионных взаимодействий и плоских групп записывается в виде 

Е Ф = У^О^пир-б)) 4 , (7) 

ф 

где ш - кратность торсионного барьера, 8 - сдвиг фазы, константы К ф - определяют 
высоты потенциальных барьеров двухгранных углов ф . 

Ван - дер - Ваальсовые взаимодействия атомов: 


Е 


и - 



( 8 ) 


где А у и В у - параметры взаимодействия; 




( * * V 2 

е і е І 
— + 


2 2 


V 


8 і е П 


В у = 2 


( * * л 6 

е і е ] 

— + — 

2 2 


л/ 8 ' 8 ! ’ 


2 * половина минимальной энергии разделения двух атомов. 


е і - глубина взаимодействия для і-го атома, 

^ ' глубина взаимодействия для |-го атома, 

% расстояние, на котором происходит взаимодействие. 

В работе исследование термодинамических параметров системы фуллерен - мо 
лекулы водорода проводилось по следующим формулам: 

Температура системы определялась через среднюю кинетическую энергию, при 
холящуюся на одну степень свободы: 


к Б Т = 


2К 

Зп - 3 ’ 


(9) 


где к Б - константа Больцмана, а К - кинетическая энергия системы фуллерен - молеку- 
лы водорода: 



( 10 ) 
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Давление вычислялось по следующей формуле: 


Р = 


1 

зж 


2 к-ХЦ-ОС/ 

и 

/</ 


(И) 


Здесь ] означает, что наряду с і -й частицей рассматриваются также все ее образы в 
соседних ячейках и выбираются координаты той частицы, которая оказалась ближе 

всего к 1 -й частице; Ж - сила, действующая на і-ю частицу со стороны ] -й; IV- объем 
расчетной ячейки [15]. 

Решение системы уравнений ( 1 ) осуществлялось численно по схеме Верлета, опи- 
сание которой представлено в [1 1,12,14,15]. 


РЕЗУЛЬТАТЫ РАСЧЕТОВ 


Рассмотрим основные закономерности поглощения и выделения водорода фуллерена- 
ми, выявленные при численном моделировании данных процессов. 

На рис. 2 приведен фуллерен С 240 в окружении молекулярного водорода при тем- 
пературах аОО (а) и 50 (б) градусов Кельвина (К 0 ). Видно, что в первом случае фулле- 
рен находится в “облаке” молекулярного водорода, а во втором - покрыт плотным сло- 
ем молекул водорода (рис. 26). 



а) б) 

Рис. 2. Фуллерен С 2 4о (1) в окружении молекулярного водорода (2) 
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Рис. 3. Зависимость количества адсорбированных молекул водорода N от темпе- 
ратуры Т для фуллерена С 180 

На рис.З представлен график зависимости количества адсорбированных молекул 
водорода N при различных температурах Т для фуллерена С 180. При моделировании 
использовалась ячейка с длиной ребра 50А с периодическими граничными условиями. 
Радиус фуллерена С 180 составляет 6А. Время процесса накопления водорода составля- 
ло 500пс. Из расчетов следует, что процесс адсорбции существенно зависит от темпе- 
ратуры. При уменьшении температуры количество молекул водорода, адсорбирован- 
ных фуллереном С 180, интенсивно возрастает. В диапазоне от 50 до 100 градусов 
Кельвина количество данных молекул возрастает в 5 раз 

Структура адсорбированного водородного облака в виде гистограмм распреде- 
ления количества молекул водорода и плотности водородного облака вокруг фуллерена 
С 180 при Т=40К представлена на рис.4, и рис.5, соответственно. Из данных гисто- 
грамм видно, что молекулы водорода концентрируются в узкий внешний слой толщи- 
ной ~5 ангстрем. Минимальное расстояние между атомами углерода молекулы фулле- 
рена и атомами молекулы водорода равно 2,4А. 

Работу фуллерена С 180 в качестве элементарной ячейки водородного аккумуля- 
тора, позволяющей накапливать и выделять водород, иллюстрирует рис. 6, на котором 
приведен график изменения содержания водорода на фуллерене С 180 при переменных 
температурах. В течение первых 220пс температура равнялась 40К, затем увеличилась 
до 300К. На первом этапе происходит накопление и “стабильное” хранение водорода 
на фуллерене. На втором этапе (при нагреве) водород отсоединяется от фуллерена и 
может быть использован в качестве источника энергии. 

В заключение анализа проведенных численных расчетов рассмотрим структуру 
водородного облака, образованного молекулами водорода вокруг фуллерита, состояще- 
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Рис. 4. Гистограмма радиального распределения молекул водорода N в зависимо- 
сти от расстояния 8 (А) от центра фуллерена С 180 при Т=40К 



Рис. 5. Гистограмма радиального распределения плотности молекул водорода 
р,кг/м 3 в зависимости от расстояния 8, А от центра фуллерена С 180 при 
Т=40К 
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Рис. 6. График адсорбции/десорбции молекул водорода фуллереном С 180 при 
переменных температурах 



Рис. 7. Структура водородного облака, образованная молекулами водорода вокруг 
фуллерита, состоящего из двух фуллеренов С 180 


го из двух фуллеренов С180 (рис.7). Как видно из данного рисунка, молекулы водорода 
неравномерно распределены вокруг фуллерита. Плотность водородного облака между 
фуллеренами намного выше его плотности в периферийной части фуллерита. 
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ВЫВОДЫ 

Проведенные расчеты позволили выявить структуру водородных образований на фул- 
лерене и в фуллеритах, и установить основные закономерности поглощения, хранения 
и выделения водорода указными наноструктурами. В целом численные исследования 
показывают, что фуллерены и фуллериты могут являться эффективными элементарны- 
ми ячейками аккумуляторов водорода, работающих в циклическом режиме: поглоще- 
ние, хранение и выделение водорода. 

Работа поддержана комплексной программой фундаментальных исследований 
Президиума РАН № 26 «Водородная энергетика», гос. контракт № 1 0002-25 1/П- 
26/117-383/290404-138. 
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ОезспЬес! аЬоѵе іурез оР Риііегепез із регРогтеО іоо. Кіпеііс оР ассишиіаііоп ргосезз оР 
Ьус1го§еп шоіесиіез із саісиіаіесі. ТЬегшосІупатіс рагатеіегз Рог зіаЫе Ьу0го§еп зіога§е 
Ьу іЬе Риііегепе аге езіаЫізЬеё. 8ішсіигез сгеаіесі Ргот Ьусіго ё еп тоіесиіез агоипсі 
Риііегепез 0игт§ §іога§е ргосезз аге Ьиііі. ТЬегшосІупатіс сопсііііопз апсі кіпеііс оР Ьу- 
ёгодеп сіезогрііоп Ргот Риііегепез аге саісиіаіесі. 
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ОЦЕНКА ГИДРАВЛИЧЕСКИХ ПОТЕРЬ В ГАЗОХОДЕ ПРИ 
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В. М. АТАМАНЮК 


Военная академия РВСН имени Петра Великого, Москва, Россия 


АННОТАЦИЯ. Для стенда сжигания при пониженном давлении твердотопливных 
зарядов снимаемых с вооружения ракет проведены расчеты для определения гид- 
равлических потерь в газоходе при впрыске охладителя и раствора нейтрализатора. 
Показано влияние скорости впрыска, расхода и угла наклона струи охладителя на 
гидравлические потери. 


Одним из возможных путей решения проблемы ликвидации твердотопливных зарядов 
снимаемых с вооружения ракет является создание стендов для их сжигания при пони- 
женном давлении (без соплового блока). Качественное проектирование и строительство 
такого стенда невозможно без соответствующего научно-технического сопровождения, 
обеспечивающего оптимизацию конструкции и основных технологических параметров 
стендового оборудования. Имеющийся к настоящему времени опыт сжигания твердо- 
топливных зарядов не позволяет оптимизировать технологические параметры истече- 
ния горячих турбулентных струй продуктов сгорания из камеры двигателя в газоход с 
целью обеспечения высокого уровня технической и экологической безопасности. В 
связи с этим представлялось целесообразным в качестве одного из этапов предвари- 
тельных работ оценить гидравлические потери в газоходе. 

Элементами, создающими гидравлические потери в газоходе, являются: 

- расширение струи от отверстия в корпусе РДТТ до канала газохода диаметром 4,5 м; 
-впрыск охладителя через отверстия на периферии стенки; введение дополнительной 
массы со своим расходом и своей ориентацией относительно направления течения га- 
за вызывает потери импульса на ускорение этой массы и, соответственно, потери дав- 
ления; 


впрыск водного раствора нейтрализатора, который действует на потери давления ана- 
логично впрыску охладителя. 


Расчет потерь при внезапном расширении 
газоходе был выполнен по аналогии с примером 


струи и новых параметров потока в 
Г. Н. Абрамовича по расчету потерь 
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давления на мерном сопле в дозвуковом трубном течении [1]. В соответствии с этим 
приведенная скорость Л гх в газоходе находится из уравнения: 


40=40 


+ 


к + 1 


к - 1 


у(А) 


кк 

Г 


где 2 (Л гх ), 2 (Л С ), у(Л с ) - газодинамические функции соответствующих приведенных ско- 
ростей. 

После подстановки значений газодинамических функций приведенной скорости 
Л с -0,604, к-1,33, г(Л)=Л+1/Л = 1,6556 и у(Л)=\,0043 (табличное значение) и отношения 
площадей поперечных сечений газохода и сопла Д г Лф= 15,9/0,4864 =32,7 получили 
?(Л г *)-52,38. Значение 4 гх =0,02 было получено в результате решения квадратного урав- 
нения газодинамической функции. 

Относительные потери полного давления при внезапном расширении дозвуково- 
го потока были найдены по зависимости 


а = К = К 40 
К р„ 40 ' 

где с/ (Л) — табличная газодинамическая функция. 

Подстановкой в эту зависимость значений функции для приведенных скоростей в 
сопле и газоходе получили а=0,7824. Если полное давление на срезе Р с =0,128 МПа, то 
после зоны рециркуляции газоходе оно составит Р гх і = 0, 128.0,7824=0, 100 15 МПа а 
средняя скорость после зоны рециркуляции ѵг«20 м/с. 

Длина зоны рециркуляции при внезапном расширении турбулентной струи приня- 
та, как в случае эжектора, равной длине, при которой расширяющаяся струя касается 
стенок газохода. Если угол расширения струи диаметром 0,787 м составляет 28-30°, а 
диаметр ограничивающего канала равен 4,5 м, то струя будет касаться стенок на рас- 
стоянии ~7 м. Разрежение в зоне рециркуляции будет сокращать ее длину. 

Расчет потерь давления при подводе массы охладителя проводили с учетом от- 
носительного расхода охладителя т * ох = т ох /т г , а также ориентации и скорости струй 
охладителя, характеризуемых параметром 

к=(\ѵ ох .со8а)/\ѵ г< 

где а- угол между направлением впрыска охладителя и направлением потока газа в га- 
зоходе. 

Относительные потери давления находили по зависимости [2] 
<у см =2т * ох (2 +т * ох -к), полученной из рассмотрения процесса втекания вторичной среды 


320 


ХИМИЧЕСКАЯ ФИЗИКА И МЕЗОСКОПИЯ. Том 7, № 3 


ОЦЕНКА ГИДРАВЛИЧЕСКИХ ПОТЕРЬ В ГАЗОХОДЕ ПРИ СЖИГАНИИ ТВЕРДОТОПЛИВНЫХ 

ЗАРЯДОВ 


в канал, элемент которого представлен на рис Л. Сносящий поток в канале до смешения 
с втекающими струями движется со скоростью \ѵ г . В сечении 1-1 в канал вводятся 
струи вторичной среды, имеющие горизонтальную составляющую скорости 
^ ох~^'охСояа. Вертикальная составляющая скорости втекающей струи, т.е. величина, 
пропорциональная расходу через отверстия т ох будет зависеть от перепада давлений на 
форсунках, при котором происходит истечение. В некотором сечении II скорости сред 
выравняются и приобретут среднее значение \ѵ г2 . 

Для выделенного элемента канала (рис.1) уравнение количества движения запи- 
шется в виде 

Р ІРг 1 “ Рг 2 ) = т г2™г2 ~ 

Следует заметить, что: 

т г 2 = т г і+т ох =т г і(1+ т* ох ); 
м> г і(1+ т* ох ); 

ГПг1= ^гіРгР 

После подстановки этих обозначений в уравнение количества движения и элемен- 
тарных преобразований было получено выражение для потерь давления на выделенном 
элементе канала в виде: 

д Ді -2 = Р 2 1 2т * ох (2 + т' ох -/(:)] 
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или в виде коэффициента относительных потерь статического давления а ш =Лр\. 2 /(р 
м> г і /2) на рассматриваемом участке в долях скоростного напора сносящего потока. Ко- 
эффициент потерь статического давления может иметь как положительные, так и отри- 
цательные значения в зависимости от соотношения величин к и т* ох . Если к= 0, т.е. вте- 
кающие струи не имеют горизонтальной составляющей скорости, то на разгон втекаю- 
щих струй до скорости м> г2 сносящий поток затратит часть своей энергии и статическое 
давление газа в сечении 2-2 будет меньшим, чем в сечении 1. В этом случае при любых 
соотношениях расходов охладителя и газа т* ох коэффициент потерь статического дав- 
ления будет иметь положительное значение, быстро возрастающее с увеличением от- 
носительного расхода охладителя (пример из [2] на рис.2). Если горизонтальная со- 
ставляющая скорости втекающей среды \ѵ' ох превышает скорость сносящего потока \ѵ гІ 
более чем в два раза, напр. к= 2,5, то до значения т* о *=0,5 коэффициент потерь статиче- 
ского давления будет иметь отрицательное значение, при т* о *=0,5 этот коэффициент 
<7см=0 и при т* ох >0,5 относительные потери давления примут положительные значения. 

При втекании струй в сносящий поток будут иметь место потери полного напора 
на смешение, что следует из уравнения Бернулли: 



С учетом того, что \ѵ г2 = \ѵ г/ (1+ т* ох ) и 



-ОД 


Рис.2. Влияние к и п^т ох на коэффициент потерь статического давления а см [2] 
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Р\~Рг = Р ™ 2 ' [ 2 »С ( 2 + "С -*)], 
после соответствующих подстановок получим: 

А Р* = Р 2 ' [ ( т і} + 2 т 'ох~ 2 ™'ох к ]• 

Семейство кривых сг * с м имеет качественно тот же характер, что и семейство кри- 
вых коэффициента потерь статического давления. 

При значениях к от нуля до единицы и при всех значениях т* ох полный напор 
сносящего потока падает, так как этот поток затрачивает энергию на разгон струи. При 
значениях к> 1 и малых значениях т* ох имеются области, где величина полного напора 
увеличивается (а* см <0) за счет разгона сносящего потока втекающими струями; однако 
при дальнейшем увеличении т* ах значение а* см переходит через ноль и становится 
снова положительным. Здесь потери энергии полного напора в сносящем потоке на по- 
ворот струй и перемешивание превосходит увеличение энергии за счет горизонтальных 
составляющих втекающих струй, разгоняющих сносящий поток. 

Эти результаты получены в предположении, что компактные струи охлаждаю- 
щей воды при поступлении в канал превращаются в пар той же самой плотности, что и 
газ. Для рассмотренной задачи в интегральной постановке не имело значения число 
форсунок, их диаметр и перепад давлений на них. Также не учитывались потери дав- 
ления на обтекание газом решетки из струй в поперечном сечении канала газохода. 

Результаты расчетов приведены в табл.1 для возможных случаев ввода воды с 
различными скоростями и под различными углами к оси канала на относительные по- 
тери статического давления. 

Принятые модули (отношения расходов воды к расходу газа) ввода охладителя на 
уровне 12-15 отражают опыт работы в основном на модельных конструкциях, причем 
без анализа к.п.д. использования воды. В этих экспериментах выброс капельно-жидкой 
воды из модельных газоходов не учитывался. 

Если известна термодинамическая энтальпия продуктов сгорания сжигаемого 
топлива, то за модуль в качестве исходной величины может быть принято значение 

т ох = Нтг/Н исп 

где Н исп - теплота испарения воды 2,26 МДж/кг при 100°С; 

#7т- энтальпия продуктов сгорания. 

Если #тт находится на уровне 25 МДж/кг, то в первом приближении 
т 0 л = 25:2,26=1 1,06 кг воды на 1 кг газа. 
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Таблица 1 . Относительные потери давления, вызванные вводом массы охладителя 


Скорость 

впрыска, 

м/с 

т* ох 

ст см при угле отклонения потока от оси газохода, град 

СГ*см 

при а=90° 

15 

30 

45 

60 

75 

90 

100 

2 

-3,316 

-1,32 

1,858 

6,0 

10,82 

16 

8 


4 

9,368 

13,36 

19,72 

28 

37,65 

48 

24 


8 

82,736 

90,72 

103,44 

120 

139,3 

160 

80 


10 

143,42 

153,4 

169,3 

190 

214,1 

240 

120 

50 

2 

6,341 

7,34 

8,929 

11,0 

13,41 

16 

8 


4 

28,682 

30,68 

33,86 

38,0 

42,82 

48 

24 


6 

67,023 

70,02 

74,79 

81,0 

88,24 

96 

48 


8 

121,36 

125,36 

131,72 

140,0 

149,65 

160 

80 


10 

151,71 

196,7 

204,65 

215,0 

227,1 

240 

120 

10 

2 

14,068 

14,268 

14,58 

15,0 

15,48 

16 

8 


4 

44,14 

44,54 

45,17 

46,0 

46,96 

48 

24 


6 

90,20 

90,80 

91,76 

93,0 

94,45 

96 

48 


8 

152,27 

153,1 

154,34 

156,0 

157,93 

160 

80 


10 

230,34 

231,3 

232,93 

235,0 

237,41 

240 

120 


Если учесть удельную теплоемкость капельно-жидкой воды при нагреве ее от 
20°С до расчетной температуры охлаждения газа, то 


т о* = тт С Л 

Н »сп+с р (Т рип -Т 0 ) 

Учитывая, что удельная теплоемкость воды составляет около 4 кДж/кг-К, величи- 
на поправки может составить около 10%. Введение в этот насыщенный водяной пар с 
небольшой примесью газа дополнительно водного раствора нейтрализатора приведет к 
выпадению капель воды. 

Ввод охлаждающей жидкости (воды) в газоход с большим поперечным сечением 
предпочтительно осуществлять компактными струями достаточно большого диаметра с 
тем, чтобы потери массы струй из-за испарения и разбрызгивания вследствие гидроди- 
намической неустойчивости были равномерными по достаточно большой длине струи, 
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соизмеримой с радиусом газохода. 

Можно предположить, что в высокотемпературном газовом потоке теплообмен 
струи с газом будет осуществляться вынужденной конвекцией в режиме пленочного 
охлаждения (нагрева), когда тепло- и массообмен происходят через паровую пленку. В 
этом случае критический удельный тепловой поток (первая критическая плотность теп- 
лового потока) для воды составляет 0, 7- 1,2 МДж/м 2 с [3]. Если поглощаемый тепловой 
поток расходуется только на парообразование с образованием потока массы с поверх- 
ности струи, то изменение по времени поперечного сечения элемента струи, движуще- 
гося со скоростью \ѵ ох , можно описать уравнением 




Так как Ѵ=лК 2 у, сІѴ = 2лК(1)сіКу и 8(1) = 2лК(і)у, то 


сШ 

(ІІ 


Ч крит 1 п _ Ч крит 1 

_ , , М/С ИЛИ Л 0 - К = і 

р я„ Р ■ 


Поскольку струя движется с постоянной скоростью, то изменение радиуса струи 
по длине можно выразить 


К = К 0 


Я крит 1 

X 

Н исп Р™ 


М. 


Дальнобойность такой струи (без учета изгиба сносящим потоком и разбрызгива- 
ния) будет равна 

/ _ К 0 Н испР™о Х 
1 ~ ' , ~ , м 

У крит 

Если начальный радиус К 0 = 1 см, теплота испарения воды 2,26 МДж/кг, плот- 
ность воды 1000 кг/м 3 , скорость струи 100 м/с, то длина струи составит 2,26 м, т.е. 
практически достигнет оси канала газохода. При принятом в проекте перепаде давле- 
ний на форсунках 0,5МПа начальная скорость струи в лучшем случае составит около 
30 м/с и дальность действия струй диаметром 20 мм ограничится пристеночной зоной 
толщиной около 70 см. Насколько известно, такие крупные форсунки в проекте не пре- 
дусмотрены, а при использовании форсунок диаметром 10 мм толщина пристеночной 
зоны ограничится 35 см. 
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В любом случае, за зоной впрыска охладителя должна располагаться буферная 
зона для выравнивания состава парогаза. Лучшие условия для перемешивания будут в 
случае использования дальнобойных и наклонных струй, чем струй, направленных по 
нормали к оси канала газохода. 

В зоне ввода нейтрализатора форсунки ввода могут быть расположены по попе- 
речному сечению канала газохода, а не только по периферии, поскольку температура 
парогаза уже низкая. 

Из материального баланса следует, что расчетное начальное давление воды на 
орошение на входе в газоход составляет 0,5 МПа, а конечное - 0,2 МПа. Это означает, 
что расход воды на орошение в течение 150с (за время сгорания заряда) падает в 


0,5 

Іо, 2 


= 2,5 = 1,58 


раза, т.е. вместо заданного значения 2000 кг/с в конце прожига в 


газоход будет поступать 1266 кг/с. Таким образом, заведомо создается высокая вероят- 
ность проскока необработанных продуктов сгорания в сепаратор. Этот проскок будет 
возрастать с увеличением времени прожига в силу каких-либо причин, например, 
вследствие снижения давления в камере сгорания. Этот недостаток обусловлен приня- 
той вытеснительной системой подачи воды для орошения с использованием воздушно- 
го аккумулятора давления без стабилизатора расхода воздуха. Расчеты параметров газа 

в емкости при ее опорожнении в различных случаях, в т.ч. и в динамике, можно найти в 
книге [4]. 

Из материального баланса также следует, что камера сгорания работает при дав- 
лении 1 МПа, что повышает скорость горения по сравнению с атмосферными условия- 
ми. При этом поджатие камеры осуществляется за счет гидравлических сопротивлений 
газохода, которые будут изменяться в зависимости от модуля т* ох . Снижение модуля к 
концу прожига приведет к снижению гидравлических потерь, снижению давления в 
камере, уменьшению скорости горения и увеличению времени прожига, т.е. опять к 
снижению модуля. 

Применение давления в газоходе выше атмосферного давления будет вызывать 
утечки необработанных продуктов сгорания через любые щели. В идеале, в газоходе 
должно поддерживаться давление ниже атмосферного. 

Принятая концентрация вводимого раствора нейтрализатора (10%-ного раствора 
каустической соды) представляется недостаточно оправданной. Это значение является 
общепринятым для нейтрализации газов, в которых вода практически отсутствует и на 
поверхности капель возможно образование слоя продуктов нейтрализации. В данном 
же случае парогаз имеет более чем 7-кратное содержание воды виде насыщенного пара 
и мелкодисперсных капель (тумана). Поэтому концентрация каустической соды может 
быть увеличена по меньшей мере вдвое и за счет этого может быть сокращено потреб- 
ление воды до 10%. При поступлении капель нейтрализатора в объем парогаза они бу- 
дут насыщаться влагой за счет диффузии и концентрация будет убывать. Ухудшение 
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ЗАРЯДОВ 


условий распыления форсунками при увеличении концентрации щелочи вряд ли будет 
наблюдаться, поскольку вязкость раствора при 20°С с 10% масс. ИаОН составляет 1,86 
сП (как вода при 0°С), а с 20% - 4,48 сП, что сравнимо с вязкостью анилина при этой 
температуре [5]. 

СПИСОК ЛИТЕРАТУРЫ 

1. Абрамович Г.Н. Прикладная газовая динамика. -М: Наука, 1976. 888 с. 

2. Михайлов А. И., Горбунов Г.М., Борисов В. В. и др. Рабочий процесс и расчет камер 

сгорания газотурбинных двигателей. -М.: Оборонгиз, 1959. 285 с. 

3. Кутателадзе С. С. Основы теории теплообмена. М-Л.: Машгиз, 1957. 383 с. 

4. Беляев Н.М., Белик Н.П., Уваров Е.И. Реактивные системы управления космических 

летательных аппаратов. -М.: Машиностроение, 1979. 232 с. 

5. Краткий справочник физико-химических величин. Под ред. Мищенко К.П., Равде- 
ля А. А. -Л.: Химия, 1967. 182 с. 


81ЛММАКУ. Ап езіітаііоп оГ гезізіапсе Ьеаб Гог Іезі ЬепсЬ оПоѵѵ ргеззиге сотЬизііоп 
оГ зоііё ргореііапі сЬаг§е оГ аппіЫІаЫе госкеіз із скте. ТЬе іпВиепсе оГ іпіеі ѵеіосііу, 
сопзитрііоп оГ \ѵаіег апсі ]е! ІіИ ап§1е оп гезізіапсе Ьеаб із зЬо\ѵп. 


ХИМИЧЕСКАЯ ФИЗИКА И МЕЗОСКОПИЯ. Том 7, № 3 


327 


РЕФЕРАТЫ 


УДК539.3:621. 454.3 

Напряженно-деформированное состояние зарядов РДТТ в области донного торца. 

Дунаев В.А., Каширкин А.А., Евланова О.А. ХИМИЧЕСКАЯ ФИЗИКА И 
МЕЗОСКОПИЯ. 2005. Т. 7. №3. С. 237-242. 

Разработана комплексная математическая модель, содержащая модели одномер- 
ного квазистационарного течения, двумерного течения газов в области торцов заряда и 
напряженно-деформированное состояние заряда. На основе данной модели проведен 
анализ движения потоков продуктов сгорания в области у переднего днища и напря- 
женно-деформированного состояния донного торца заряда. Ил. 2. Библиогр. 5. 

УДК 621.454.3 

Баллистические характеристики активно-реактивных снарядов с применением 
композиционных материалов и надежность их маршевых ракетных двигателей. 

Александров Д.П., Захаренков В.Ф. ХИМИЧЕСКАЯ ФИЗИКА И МЕЗОСКОПИЯ 2005 
Т. 7. №3. С. 243-250. 

Изложен машинный эксперимент по изучению влияния механических характе- 
ристик твердого топлива и материалов корпусов маршевых ракетных двигателей на 
твердом топливе (РДТТ) при замене металлических корпусов на стеклопластиковые и 
боралюмоуглепластиковые на боевую эффективность активно-реактивных снарядов 
(АРС) и надежность их функционирования. Даны статистический анализ работоспо- 
собности РДТТ и рекомендации по повышению технических характеристик изделий. 
Табл. 2. Ил. 4. Библиогр. 4. 

УДК 539.3:623.4 

Исследование взаимодействия комбинированных ударников с преградами. 

Глазырин В.П., Орлов Ю.Н., Орлов М.Ю. ХИМИЧЕСКАЯ ФИЗИКА И 
МЕЗОСКОПИЯ. 2005. Т. 7. №3. С. 251-258. 

Разработана методика компьютерного моделирования взаимодействия ударни- 
ков с преградами из различных конструкционных материалов. Спрогнозированы ре- 
зультаты взаимодействия, а также выявлено влияние материалов сердечников на про- 
бивное действие ударников. Установлено, что за счет изменения компоновки и формы 
ударника для дозвуковой скорости взаимодействия можно существенно увеличить его 
пробивное действие. Табл. 3. Ил. 4. Библиогр. 5. 


ХИМИЧЕСКАЯ ФИЗИКА И МЕЗОСКОПИЯ. Том 7, № 3 


УДК 677-486.1 

Особенности поведения углерод-углеродных материалов при ударно- 
волновом нагружении. Апакидзе Ю.В., Тихонов А.А., Уткин А. В., Гафаров Б.Р. 
ХИМИЧЕСКАЯ ФИЗИКА И МЕЗОСКОПИЯ. 2005. Т. 7. №3. С. 259-265. 

Приведены результаты экспериментальных исследований поведений углерод- 
углеродных материалов при ударно-волновом нагружении. Полученные результаты по- 
казывают высокую чувствительность структуры фронта ударных волн и откольной 
прочности к структуре композиционных материалов. Ил. 4. Библиогр. 8. 

УДК 623.4.017 

Термоэрозионная и механическая стойкость стволов артиллерийских автоматов. 

Александров Д.П., Захаренков В.Ф., Цепелев В.С. ХИМИЧЕСКАЯ ФИЗИКА И 
МЕЗОСКОПИЯ. 2005. Т. 7. №3. С. 266-276. 

Приводится математическое описание модели термоэррозионной и механиче- 
ской стойкости стволов артиллерийских автоматов в условиях стрельбы короткими или 
длинными очередями с темпом стрельбы 600 выстр./мин. Комплексная математическая 
модель включает описание процессов внутренней баллистики, теплообмена, теплопро- 
водности и износа с учетом структурно-фазовых превращений в материале орудийной 
стали, а также механических и температурных напряжений 1 -го и 2-го рода, возникаю- 
щих при циклическом высокоимпульсном тепловом нагружении. Результаты числен- 
ных исследований по разработанной модели иллюстрируются примером расчетов 30- 
мм гипотетической автоматической пушки, имеющей начальную скорость снаряда 960 
м/с и максимальное давление 380 МПа. Ствол изготовлен из орудийной стали 
38ХНМФА. Ил. 10. Библиогр. 10. 

УДК 539.2: 539.3 

Методика расчета упругих параметров наноэлементов. Вахрушев А.В.. 
Шушков А. А. ХИМИЧЕСКАЯ ФИЗИКА И МЕЗОСКОПИЯ. 2005. Т.7. №3. С. 277-285" 

Разработана методика расчета упругих констант нано элементов методом согла- 
сования перемещений составляющих их атомов с полями перемещений упругих 
эквивалентных элементов при различных видах нагружения. Приведены результаты 
расчетов модуля упругости наночастиц, нагруженных растягивающими сосредото- 
ченными силами. Получены расчетные зависимости модуля упругости Юнга от диа- 
метра наночастиц. Ил.7. Библиогр. 14. 


ХИМИЧЕСКАЯ ФИЗИКА И МЕЗОСКОПИЯ. Том 7, № 3 


329 


УДК 532.519.6 

Численное исследование течений жидкости с переменными вязкими свойствами в 
изогнутых каналах. Вахрушев А. А., Липанов А.М., Вахрушев А. В. ХИМИЧЕСКАЯ 
ФИЗИКА И МЕЗОСКОПИЯ. 2005. Т.7. №3. С. 286-300. 

Рассмотрена задача моделирования двумерных течений неньютоновских жидко- 
стей в криволинейных каналах. Проведены численные расчеты параметров псевдопла- 
стической среды, описываемой с помощью реологической модели Пауэлла - Эйринга, и 
исследованы физические явления, возникающие при локальном изменении температу- 
ры стенки канала. Ил. 6. Библиогр. 8. 

УДК 539.3:621.454.3.004.82 

О критерии безопасного сжигания дефектных зарядов твердого топлива при 
пониженных давлениях. Поспелов Д.А., Мелешко В.Ю., Наумов П.Н., Павловец Г Я 
ХИМИЧЕСКАЯ ФИЗИКА И МЕЗОСКОПИЯ. 2005. Т.7. №3. С. 301-308. 

Рассмотрены дефекты твердотопливных зарядов снимаемых с вооружения ра- 
кет, влияющие на процесс сжигания. Предложены критерии, позволяющие оценить 
влияние дефектов заряда на увеличение поверхности горения и расхода продуктов сго- 
рания. Определены условия, при которых растрескивание внутренней поверхности за- 
ряда можно считать не опасным для его утилизации. Библиогр. 3. 

УДК 539.2: 544.723 

Моделирования процессов поглощения и выделения водорода фуллеренами и 
фуллеритами. Липанов А.М., Вахрушев А.В., М.В.Суетин, А. А. Вахрушев 
ХИМИЧЕСКАЯ ФИЗИКА И МЕЗОСКОПИЯ. 2005. Т.7. № 3. С. 309-318. 

Методом молекулярной динамики выполнены расчеты адсорбции, хранения и 
выделения водорода фуллеренами С 20 , С 60 , С 80 , С 180> С 240 и С 540 и сформированными из 
фуллеритами при различных внешних условиях. Рассчитаны кинетика адсорбции мо- 
лекул водорода фуллеренами; определены термодинамические параметры, необходи- 
мые для устойчивого хранения водорода фуллереном; построена структура, образован- 
ная молекулами водорода на фуллерене при хранении; определены термодинамиче- 
ские условия выделения водорода. Ил.7. Библиогр. 15. 


330 


ХИМИЧЕСКАЯ ФИЗИКА И МЕЗОСКОПИЯ. Том 7, № 3 


УДК 532.55:621.454.3.004.82 

Оценка гидравлических потерь в газоходе при сжигании твердотопливных 
зарядов. Карелин В. А., Мелешко В.Ю., Шустачинский В. С., Атаманюк В.М. 
ХИМИЧЕСКАЯ ФИЗИКА И МЕЗОСКОПИЯ. 2005. Т.7. №3. С. 319-327. 

Для стенда сжигания при пониженном давлении твердотопливных зарядов сни- 
маемых с вооружения ракет проведены расчеты для определения гидравлических по- 
терь в газоходе при впрыске охладителя и раствора нейтрализатора. Показано влияние 
скорости впрыска, расхода и угла наклона струи охладителя на гидравлические потери. 
Табл. 1. Ил. 2. Библтогр. 5. 
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